arXiv:1506.05297v3 [cs.SY] 22 Aug 2016 


DECENTRALIZED ABSTRACTIONS FOR MULTI-AGENT SYSTEMS 
UNDER COUPLED CONSTRAINTS 

D. BOSKOS AND D. V. DIMAROGONAS 


Abstract. The goal of this report is to define abstractions for multi-agent systems with 
feedback interconnection in their dynamics. In the proposed decentralized framework, we 
specify a finite or countable transition system for each agent which only takes into account 
the discrete positions of its neighbors. The dynamics of the considered systems consist 
of two components. An appropriate feedback law which guarantees that certain system 
and network requirements are fulfilled and induces coupled constraints, and additional free 
inputs which we exploit in order to accomplish high level tasks. In this work, we provide 
sufficient conditions on the space and time discretization for the abstraction of the system’s 
behaviour which ensure that we can extract a well posed and hence meaningful transition 
system. Furthermore, these conditions include design parameters whose tuning provides the 
possibility for multiple transitions, and hence, enable the construction of transition systems 
with motion planning capabilities. 


1. Introduction 

High level task planning for multi-agent systems constitutes an active area of research which 
lies in the interface between computer science and modern control theory. A challenge in this 
new interdisciplinary direction constitutes of the problem of defining appropriate abstractions 
for continuous time multi-agent control systems, which can be used as a tool for the analysis 
and control of large scale systems and the synthesis of high level plans HU, |2U. Robot motion 
planning and control constitutes a central field where this line of work is applied llOj . [5,. In 
particular, the use of a suitable discrete system’s model allows the automatic synthesis of dis¬ 
crete plans that guarantee satisfaction of the high level specifications. Then, under appropriate 
relations between the continuous system and its discrete analogue, these plans can be converted 
to low level primitives such as sequences of feedback controllers, and hence, enable the contin¬ 
uous system to implement the corresponding tasks. Such tasks in the case of multiple mobile 
robots in an industrial workspace could include for example the following scenario. Robot 1 
should periodically visit regions A , B, while avoiding C and after collecting an item of type X 
from robot 2 at location D, store it at location E. 

In order to synthesize high level plans, it is required to specify an abstraction of the original 
system, namely a system that preserves some properties of interest of the initial system, while 
ignoring detail. Results in this direction for the nonlinear single plant case have been obtained 
in the papers |23| and |32| . where the notions of approximate bisimulation and simulation are 
exploited for certain classes of nonlinear systems, under appropriate stability, and completeness 
assumptions, respectively. The notion of bisimulation, has its origin in computer science jjj, 
and guarantees that if the initial system and its abstraction are bisimilar, then the task of 
checking feasibility of high level plans for the original system reduces to the same task for its 
abstraction and vice versa. Bisimulation relations between transition system models of discrete 
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or continuous time linear control systems with finite affine observation maps were explicitly 
characterized and constructed in |[22], providing also a generalization of the notion of state 
space equivalence between continuous time systems [50] . 

Another abstraction tool for a general class of systems is the hybridization approach [5], where 
the behaviour of a nonlinear system is captured by means of a piecewise affine hybrid system on 
simplices. Motion planing techniques for the latter case have been developed in the recent works 
El, El, which are also based on the abstraction and controller synthesis framework provided 
in El, El, and further studied in [7]. Other abstraction techniques for nonlinear systems 
include |25| , where discrete time systems are studied in a behavioral framework, the sign based 
abstraction methodology introduced in [29], which is based on Lie-algebraic type conditions 
and PP, where box abstractions are studied for polynomial and other classes of systems (for a 
literature survey on the subject see also the monograph HU). It is also noted that certain of 
the aforementioned approaches have been extended to switched ,14], [15] and networked control 
systems [31]. Furthermore, abstractions for the case of discrete time interconnected systems 
that are described by coupled difference equations, can be found in [2S] and [24], for stabilizable 
linear systems, and incrementally input-to-state stable nonlinear systems, respectively. Finally, 
we note that the control design which we adopt for the construction of the symbolic models is 
in part related to the notion of In-Block Controllability 0, [113- 

In particular, we focus on multi-agent systems and assume that the agents’ dynamics consist 
of feedback interconnection terms, which ensure that certain system properties as for instance 
connectivity or (and) invariance are preserved, and free input terms, which provide the ability 
for motion planning under coupled constraints. In this report, we generalize the corresponding 
results of our recent work [5], where sufficient conditions for well posed abstractions of the 
multi-agent system are provided for the case where the agents’ workspace is R n . A well posed 
abstraction refers to the property that the discrete state transition system which serves as 
an abstract model of the multi-agent system has at least one outgoing transition from each 
state. The extension in this work is twofold. First, the results on admissible space and time 
discretizations in [5] which ensure well posed abstractions, are now valid when the agents’ 
workspace is a general domain D of K™, provided that D is invariant for the dynamics of the 
system. Also, the corresponding framework is extended for motion planning, and sufficient 
conditions are provided which guarantee that each agent can perform multiple transitions from 
each initial discrete state. 

The rest of the report is organized as follows. Basic notation and preliminaries are introduced 
in Section 2. Section 3 is devoted to the formulation of the problem and motivates the control 
design that will be utilized for the derivation of the symbolic models. In Section 4, we define 
well posed abstractions for single integrator multi-agent systems by means of hybrid feedback 
controllers and prove that the latter provide solutions consistent with the design requirement 
on the systems’ free inputs. Section 5 is devoted to specific properties of the control laws 
that realize the transitions of the proposed discrete system’s model. In Section 6 we quantify 
space and time discretizations which guarantee well posed transitions with motion planning 
capabilities. The framework is illustrated through an example in Section 7 including simulation 
results. Finally, we conclude and indicate directions of further research in Section 8. 

2. Preliminaries and Notation 

We use the notation |cc| for the Euclidean norm of a vector x £ R". For a subset S of 
R n , we denote by cl(S), int(S) and dS its closure, interior and boundary, respectively, where 
dS := cl(S') \ int(S). Given R > 0 and x £ R n , we denote by B(R) the closed ball with 
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center 0 £ R n and radius R, namely B(R) := {x £ R ra : |x| < R} and B(x;R) := {y £ 
R” : \x — y\ < i?}. For two nonempty sets A, B C R” their Minkowski sum is given as 
A + B := {x + y £ R" : x £ A,y £ B}. Also, for a nonempty set A C R n its diameter is defined 
as diam(A) := sup{|x — y\ : x,y £ A}. Given a real number a £ R>o we denote by \a\ the 
integer part of a (|_aj := max{n £ N U {0} : n < a}). Finally, given a function / : X —► Y and 
a subset W of X, the notation f\w is used for the restriction of / to W. 

Consider a multi-agent system with N agents. For each agent i £ Af := {1,..., N} we use 
the notation A/* for the set of its neighbors and Ni := |A/) for its cardinality. We also consider 
an ordering of the agent’s neighbors which is denoted by ji, ■ ■ ■ ,jN i} and define the IVj-tuple 
j(i) = (j i,... Whenever it is clear from the context, the argument i in the latter notation 

will be omitted. Given an index set X and an agent i £ Af with neighbors ji,, j yv, £ A/”, 
define the mapping pr.j : X N -A X Ni+1 which assigns to each AZ-tuple (Zi,..., Z jv) £ X N the 
Ni + 1-tuple (li, lj 1 ,..., lj N .) £ X Ni+1 , i.e., the indices of agent i and its neighbors. 

We proceed by providing a formal definition for the notion of a transition system (see for 
instance [3], [22], [23]). 

Definition 2.1. A transition system is a tuple TS := ( Q,Act,—>), where: 

• Q is a set of states. 

• Act is a set of actions. 

• —> is a transition relation with —Q x Act x Q. 

The transition system is said to be finite, if Q and Act are finite sets. We also use the (standard) 
notation q — q' to denote an element (g, a, q') £—X For every q £ Q and a £ Act we use the 
notation Post(g;a) := {q' £ Q : ( q,a,q') £—>}. The transition system is called deterministic 
if for each q £ Q and a £ Act , q q' and q — q" implies that q' = q". 

3. Problem Formulation 

We focus on multi-agent systems with single integrator dynamics 

ii = Ui,Xi £ R",i £ Af (3.1) 

and consider as inputs decentralized control laws of the form 

Ui = fi(xi,Xj) + Vi,ieJ\f, (3.2) 

with Xj(= Xj(j)) := (xj lt ... ,Xj N .) £ R^" (see Section 2 for the Notation consisting of 

two terms: a feedback term /)(•) which depends on the states of i and its neighbors, and an 
extra input term «j, which we call free input. We assume that for each i £ Af it holds Xj £ D 
where D is a domain of R™ and that each /*(■) is locally Lipschitz. We also assume that the 
feedback terms /$(■) are globally bounded, namely, there exists a constant M > 0 such that 

\fi(xi,Xj)\ < M,V(xj,Xj) £ D Ni+1 . (3.3) 

Furthermore, we consider piecewise continuous free inputs that satisfy the bound 

\Vi(t)\ < Vmax,Vt > 0,i £ Af. (3.4) 

In the subsequent analysis, it is assumed that the maximum magnitude of the feedback terms 
is higher than that of the free inputs, namely, that 

Wiax < M. (3.5) 
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This assumption is motivated by the fact that we are primarily interested in maintaining the 
property that the feedback is designed for, and secondarily, in exploiting the free inputs in or¬ 
der to accomplish high level tasks. A class of multi-agent systems of the form (3.1)-(3.2) which 
justifies this assumption has been studied in our companion work [B], In particular, sufficient 
conditions are provided, which guarantee both connectivity of the network and forward invari¬ 
ance of the system’s trajectories inside a given bounded domain, for an appropriate selection of 
Umax in (3.4) which necessitates D max to satisfy (3.5). The latter forward invariance property is 
defined in the Invariance Assumption (IA) below, which we assume that the multi-agent system 
(3.1)-(3.2) satisfies for the rest of the report. 

(IA) For every initial condition x(0) £ D N and any piecewise continuous input v = (iq, ... , v n ) : 
R>o -» M. Nn satisfying (3.4), the (unique) solution of the system (3.1)-(3.2) is defined and re¬ 
mains in D n for all t > 0. < 


This assumption does not restrict the class of systems under consideration since it is satisfied 
for any forward complete system when D = R n . Recall that the system (3.1)-(3.2) is forward 
complete (see e.g, H) if for each initial condition in R n and each measurable locally essentially 
bounded input v = (iq,... ,v n ) : R>o —> R Wn , its solution exists for all positive times. Also, 
notice that due to the above bounds on the dynamics and the free input terms, the system 
(3.1)-(3.2) is forward complete. Finally, when D is bounded, as is the case in [5], a finite 
partition of the workspace by bounded sets can lead to a finite transition system which captures 
the properties of interest of the multi-agent system and hence enables the investigation for 
computable solutions with respect to high level specifications. 


In what follows, we consider a cell decomposition of the state space D (which can be regarded 
as a partition of D) and a time step St > 0. We will refer to this selection as a space and 
time discretization. For the definition of a cell decomposition we adopt a modification of the 
corresponding definition from El p 129-called cell covering]. 


Definition 3.1. Let D be a domain of R". A cell decomposition S = {S)}j e x of D , where 
I is a finite or countable index set, is a family of nonempty connected sets Si, l £ I, such that 
sup{diam(5;), l £ 1} < oo, cl(int(S))) = Si for all l £ I, int(S)) D int(<Sj) = 0 for all l ^ 1 and 
U iczSi = D. <1 


Given a cell decomposition S := {Si}i & x of D, we use the notation 1; = (. k,lj 1 ,..., lj N . ) £ I Ni+1 
to denote the indices of the cells where agent i and its neighbors belong at a certain time instant 
and call it the cell configuration of agent i. Similarly, we use the notation 1 = (l\,..., Ijy) £ I N 
to specify the indices of the cells where all the N agents belong at a given time instant and 
call it the cell configuration (of all agents). Thus, given a cell configuration 1, it is possible 
to determine the cell configuration of agent i as 1, = pr^ (1) (see Section 2 for the definition of 

pl(-))- 

Through the space and time discretization we aim at capturing reachability properties of the 
original continuous time system, by means of a discrete state transition system. Informally, we 
would like to consider for each agent i, its individual transition system with state set the cells 
of the state partition, actions defined to be all possible cells of its neighbors, and transition 
relation specified as follows. Given the initial cells of agent i and its neighbors, it is possible 
for i to perform a transition to a final cell, if for all states in its initial cell there exists a free 
input, such that its trajectory will reach the final cell at time St, for all possible initial states of 
its neighbors in their cells, and their corresponding free inputs. Feasibility of high level plans 
requires the corresponding system to be well posed (meaningful), which implies that for each 
initial cell it is possible to transit to (at least) one final cell. 
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We next illustrate the concept of a well posed space-time discretization, namely, a discretiza¬ 
tion which generates for each agent a meaningful transition system in accordance with the 
discussion above. Consider a cell decomposition as depicted in Fig. 1 and a time step St. The 
tips of the arrows in the figure are the endpoints of agent’s i trajectories at time St. In both 
cases in the figure we focus on agent i and consider the same cell configuration for i and its 
neighbors. However, we consider different dynamics for Cases (i) and (ii). In Case (i), we 
observe that for the three distinct initial positions in cell S /., it is possible to drive agent i to 
cell Sir at time St. We assume that this is possible for all initial conditions in this cell and 
irrespectively of the initial conditions of i’s neighbors in their cells and the inputs they choose. 
We also assume that this property holds for all possible cell configurations of i and for all the 
agents of the system. Thus we have a well posed discretization for system (i). On the other 
hand, for the same cell configuration and system (ii), we observe the following. For three dis¬ 
tinct initial conditions of i the corresponding reachable sets at St, which are enclosed in the 
dashed circles, lie in different cells. Thus, it is not possible given this cell configuration of i to 
find a cell in the decomposition which is reachable from every point in the initial cell and we 
conclude that discretization is not well posed for system (ii). 


Xi(St) 


,-r Xi{St) 


1 

ft 



to 




• x jl 


1 

Ci 






System (i): System (ii): 

Xi i.Xi,Xj 1 ,Xj 2 ) -f- L^,(i) Xi ,Xj 2 ) T 

Figure 1. Illustration of a space-time discretization which is well posed for 
system (i) but non-well posed for system (ii). 


One main challenge in the attempt to provide meaningful decentralized abstractions for 
system (3.1)-(3.2) is the interconnection between the agents through the _/)(•) terms. The 
latter in conjunction with the considerations above, motivates the design of appropriate hybrid 
feedback laws in place of the Vi’s which will guarantee the desired well posed transitions. 


We next define the particular feedback laws that are utilized for the construction of the 
symbolic models in this report, which also motivate the notion of well posed discretizations 
that will be formulated in the next section. Consider a cell decomposition S = {S)};ex of D 
and a time step St. For each agent i £ AT and cell configuration 1, = ((,, ., lj N .) £ X Ni+1 

of i let 


{xi,G,Xj,c) £ <Szi x {Si jl x • • • x Si jN ) 


(3.6) 


be an arbitrary iV,; + 1-tuple of reference points and define the family of feedback laws : 
[0,T(xio,Wi)) x D Ni+1 —> R" parameterized by xm £ and Wi £ W as 


kijtfaxi^jixio^i) : = ki : i itl (xi,Xj) 

+ kiA i}2 (Xio) + ki t l it 3(t;XiO,Wi), 


(3.7) 
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where 

W := -B(Ai> max ), A G (0,1), 

(3.8) 

and 

**,14,1 (*i> x i) : = fi(xi,Xj,a) - fi{xi,Xj), 

(3.9) 


(3.10) 

ki,u,3(t-,x i0 ,Wi) := F itU (xi{t)) +Wi- F itU (xi(t) +twi 
+ (! - JtXZiO - Zi,G))> 

(3.11) 

T(x i0 ,Wi) := sup{f G [0, T max ) : Xi(t) +tw t 

+ (i - ii) (x i0 - Xi,c) e Dyt g [o,t]}, 

(3.12) 

t G [0,T(x i0 , w^), (x it Xj) G D Ni+1 ,Xi o G Si^Wi G W. 

The function Fi^.(-) in (3.11) is defined as 

*1,1,0*) • fi ^ D 

(3.13) 

and Xi{') i n (3.111, (3.12) is the solution of the initial value problem 

Xi = = Xi,G, 

(3.14) 


which is defined and remains in D on the maximal right interval [0, T max ) (T max is the same 
as in (3.12)). In particular, Xi(') constitutes a reference trajectory, whose endpoint agent i 


should reach at time St , when the agent’s initial condition lies in S). and the feedback fcj i 4 (- 
above is applied (plus some extra hypotheses for the rest of the agents). The time T(x.io,Wi 


in (3.12) stands for the right endpoint of the maximal right interval for which a modification of 


the reference trajectory that depends on Xio and w-i is guaranteed to remain inside the domain 


D. Also, the parameter A in (3.8) stands for the portion of the free input that can be exploited 


for motion planning. In particular, each vector Wi from the set W in (3.8) provides a possible 
“constant velocity” of a motion that we superpose to the reference trajectory Xi{') of agent i, 
allowing thus the agent to reach all points inside a ball with center the position of the reference 
trajectory at time St. Note that the control laws fci,i ; (-) are decentralized, since they only use 
information of agent i’s neighbors states. In addition, they depend on the cell configuration 1,; 


through the reference points (x^GjNAg) which are involved in (3.9)-(3.11). 


We proceed by providing some extra intuition for the selection of the control laws in (3.9), 
(3.10) and (3.11), based on Fig. 2 below. Consider a cell decomposition of D , a time step 


St and select an agent i, a cell configuration of i and a tuple of reference points as in (3.6). 


The reference trajectory of i is obtained from (3.14), by “freezing” agent i’s neighbors at their 


corresponding reference points through the feedback term i,,i(-)- Also, by selecting a vector 
Wi in W and informally assuming that we can superpose to the reference trajectory the motion 
of i with constant speed Wi, namely, move along the curve ah(-) defined as 


Xiit) := Xi(t) + tWi,t G [0,T max ), 


(3.15) 


we can reach the point x inside the dashed ball at time St from the reference point x^g, as 
depicted in Fig. 2. In a similar way, it is possible to reach any point inside this ball by a 
different selection of . This ball has radius 


v :— Au max df, 


(3.16) 
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namely, the distance that the agent can cross in time St by exploiting the part of the free 
input that is available for planning. Our abstraction requirement is that the transition to each 
cell which has nonempty intersection with B(xi(St);r) is well posed. This will be verified for 
system ( 3 . 1)-( 3.2 ) through the establishment of condition ( 6.231 in Theorem 6.3 (Section 6 ) 
for appropriate space-time discretizations. For instance, in order for the transition to the cell 
where the point x belongs to be well posed, we will require the following. That the feedback law 
fcj i. (•), which is selected in place of agent’s i free input, guarantees that for each initial condition 
of i in cell S^, the endpoint of i’s trajectory will coincide with the endpoint of the curve Xi(-). 
In order to compensate for the deviation of the initial state with respect to the reference point 
and reach the point x , we use the extra terms and feei,, 3 (•)■ These enforce the agent 

to move with the velocity of the reference trajectory plus two constant velocity terms, one 
analogous to the displacement between the agents initial state and the reference point, and the 
other analogous to the distance between x and the endpoint of Xi(')- 



. Xi(t) 

— Xi(t ) 

Xi(St) -- B( Xi (St);r) 

Xi(8t) = x 

-- Xi($t)+tWi 


Figure 2. Consider any point x inside the ball with center Xi(St). Then, for 
each initial condition Xm in the cell Si t , the endpoint of agent’s i trajectory 
Xi(-) coincides with the endpoint of the curve ah(-), which is precisely x, and 
lies in Sp., namely, Xi(St) = Xi{8t) = x £ Sp. 


4. Abstractions for Multi-Agent Systems 


In this section we formalize the discussion in Section 3, by exploiting a class of hybrid feedback 


laws containing the control laws introduced in (3.71. One reason for employing the subsequent 


analysis in an abstract framework is that the selection in (3.9)-(3.11) is not the only possible. 


For instance, it can be shown that selecting the control laws (t ?: , x ? -) := fj(xi q, x j.c) 

/, i.r, . x .). ki t \ u 2(xio) as before and fci,i 4 , 3 (iOj) := Wi, can also provide well posed discretizations 


for system (3.1)-(3.2). However, the latter will necessitate finer admissible discretizations and 


hence, increase the complexity of the symbolic model. In the sequel, given a time step St and the 
bounds M and u max on the feedback and input terms provided by (3.3) and (3.4), respectively, 
it is convenient to introduce the following lengthscale 

^max • — St(M + U max ). (4.1) 

It follows from B B 
can travel within time St. 


(3.3), (3.4) and (4.1) that R max is the maximum distance an agent 


Before defining the notion of a well posed space-time discretization we define the class of 
hybrid feedback laws which are assigned to the free inputs Vi in order to obtain meaningful 
discrete transitions. For each agent, these control laws are parameterized by the agent’s initial 
conditions and a set of auxiliary parameters belonging to a nonempty subset W of R n . These 
parameters, as discussed in the previous section, are exploited for motion planning. In par¬ 
ticular, for every agent i, each vector £ W is in a one-to-one correspondence with a point 
inside a reachable ball for i, and the agent can reach this point by selecting the control law 
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corresponding to the specific parameter Wi. The latter provides the possibility for the agent to 
perform transitions to different cells, namely, all cells which have nonempty intersection with 
that ball. Furthermore, we note that in accordance to the control laws introduced in (3.7) for 
each agent i, the feedback laws in the following definition depend on the selection of the cells 
where i and its neighbors belong. One basic requirement for this class of controllers consists 
of conditions that guarantee well posed solutions for the system (condition (PI) in Definition 


4.1 below). We also impose a consistency requirement (condition (P2) in Definition 4.1) that 

when the states 


their magnitude does not exceed the maximum bound on the free inputs ( |3.4[ ) 
of the agent and its neighbors lie in an appropriate inflation of their corresponding cells 
particular, an overapproximation of their reachable states over the time interval [0, St]). 


Definition 4.1. Given a cell decomposition S = {S)}; 6 x of D , a time step St and a nonempty 
subset W of K", consider an agent i £ J\T and an initial cell configuration 1, = (l t , lj lf ..., ( 7A , ) 
of i. For each Xio £ Si t and w t £ W, let T(xio,Wi) > 0 and consider a mapping (•; Xw, Wi) : 
[0,T(a;jo,Wi)) x D Ni+1 —» R ra , parameterized by Xm € Su and Wi £ W. We say that fci,i ; (-) 
satisfies Property (P), if the following conditions are satisfied. 

(PI) For each xm £ Si i and Wi £ W, the mapping (•; ®jo, Wi) is locally Lipschitz continuous. 
(P2) It holds 

I(^> -Gi-Xy, Xio , wi )| ^ v max ,\/t £ [0,(5t] n [0, T{xiQ,Wi)), 

%i £ ( S+ B(R max )) D D,Xj m £ (S hm + B(R max )) n D, 

m= 1,..., Ni, x i0 £ Si t ,Wi £ W, (4.2) 


with Umax as given in ( |3.4| ) and R max as in (4.1). 

(P3) It holds T(xio,Wi) > St, for all rr i0 £ S^, Wi £ W. 


d 


The motivation for considering the time interval [0, T(xio, w i)) i n Definition 4.1 comes from 
the maximal right interval on which the modification of agent’s i reference trajectory in (3.12) 
remains inside the domain D. We next provide an extra Condition (C) for the feedback laws 
provided in the above definition, which is needed in order to define well posed discretizations. 


Definition 4.2. Consider a cell decomposition S — {S)}/ e x of D , a time step St and a nonempty 
subset W of K n . Given an agent i £ AT, a cell configuration 1, = (i i5 / 7l ,..., lj N . ) of i, a control 
law 

Vi = k it i i {t,Xi,-x. j -,Xio,w i ) ( 4 . 3 ) 

as in Definition | 4 . 1 | that satisfies Property (P), a vector Wi £ W, and a cell index l\ £ I, we say 
that Condition (C) is satisfied, or specifically, that 1.;, fci,p(-)> Wi, l\ satisfy Condition (C), if 
the following hold. For each initial cell configuration 1 with pr 7 (l) = lj, 1 = (Zj_,..., Z jv), an d for 
all £ £ M \ {i} and feedback laws 

ve = ke t \ t (t,xe,-x.j(ey,xto,we), ( 4 . 4 ) 


that satisfy Property (P) (with le = pr f (l)), the solution of the closed-loop system (3.1 )-( ph2[ ), 
(4.3)-(4.4) is well defined on [0, St] and satisfies Xi(St,x(0)) £ Si>, for all initial conditions 
x(0) £ D n with £i(0) = xm £ Si { xe(0) = £ Si e , £ £ J\f \ {*} and wp £ W, £ £ Af\ {?’}. < 


Notice that when Condition (C) is satisfied, agent i is driven to cell Sy. precisely in time 
St under the feedback law fcyi ; (-) corresponding to the given parameter Wi in the definition. 
In particular, Condition (C) ensures that the latter holds for any choice of feedback laws in 
place of the other agents’ free inputs, as long as these control laws satisfy Property (P). We 
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next provide the definition of a well posed space-time discretization. This definition formalizes 
our discussion on the possibility to assign a feedback law to each agent, in order to enable a 
meaningful transition from an initial to a final cell. 


Definition 4.3. Consider a cell decomposition S = {<S)}/gz of D , a time step St and a nonempty 
subset W of K". 


(i) Given an agent i £ A f, an initial cell configuration 1; = Zj, ,. .., lj N .) of i and a cell index 

I' £ I we say that the transition li -^A l\ is well posed with respect to the space-time 
discretization S — St, if there exist a feedback law Vi = (•; t,;o, w,) as in Definiti on i 

that satisfies Property (P), and a vector Wi £ W, such that Condition (C) in Definition 4.2 
fulfilled. 


4.1 


is 


(ii) We say that the space-time discretization S — St is well posed, if for each agent i £ Af 
and cell configuration 1, = (li, lj 1 ,..., lj N ) of i, there exists a cell index l( £ T such that the 

transition li —^A l\ is well posed with respect to S — St. 


Given a space-time discretization S — St and based on Definition |4.3[ i), it is now possible to 
provide an exact definition of the discrete transition system which serves as an abstract model 
for the behaviour of each agent. 

Definition 4.4. For each agent i, its individual transition system TSi := (Qi,Acti, — >i) 
is defined as follows: 

• Qi :=X (the indices of the cell decomposition) 

• Acti := I Ni+1 (the set of all cell configurations of i) 

• li — S-,; l( iff li —A l( is well posed, for each l i: /' £ Qi and lj = (li, lj,,, lj N ) £ Acti. < 


We have preferred to use the term actions instead of labels (as for instance in m) for the 
elements of the set Acti, because the cell configuration of i indicates how the feedback term 
/)(•) acts on and affects the possible transitions of agent i. 


Remark 4.5. (i) Given a well posed space-time discretization S — St and an initial cell con¬ 
figuration 1 = (Z 1; ..., Zjv), it follows from Definitions 4.3 and 4.4 that for each agent i £ Af it 
holds Postj(Zj; pr i (l)) ^ 0 (Postj(-) refers to the transition system TSi of each agent-see also 
Section 2). 

(ii) According to Definition 


4.3 


given a control law ki i.(-) it is possible to perform transitions 
to different cells by an alternative selection of w t . 

(iii) In addition, it is also possible to obtain different transitions by choosing an alternative 
control law. In particular, it is possible for the control laws considered in (3.7) to obtain a 
different reference trajectory in (3.14) by selecting another set of points (xi t G,*j,G) and hence, 
reach a ball which intersects different cells (see Fig. 2). <1 


Assume a well posed space-time discretization S — St is given. Based on Definition |4.3| 
and Remark |4.5| we proceed by providing a modification of each transition system TS, that 
captures additional information on the control actions that realize the individual transitions. 
In particular, for each i £ Af and cell configuration I, = (li, lj,,..., lj N ) of i we pick a control law 
ki : ii(-) which generates at least one well posed transition, i.e., such that 1,;, p(-), Wi, l( satisfy 

Condition (C) for certain Wi £ W and l( £ I (this is always possible since the discretization is 
well posed) and define for all l £ I 

K] (li)0 := {w £ W : li, fcj i i i ( , ) J w, l satisfy Condition(C)}. 


(4.5) 
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Based on (4.51, we next provide for each agent the controlled version of its individual transitions 
system. 


Definition 4.6. Consider a well posed space-time discretization and select for each agent i 
and cell configuration 1, a control law ^ (•) which generates at least one well posed transition. 
Then, the controlled individual transition system TSf := ( Q i: Act —>?) of each agent i 
is defined as follows 


• Qi~X 

• Acti := I Ni+1 x 2 W 




l'i and [wi] = ^ 0, for each kJl G Q u 1, 


k —►? l[ iff k — 

X Ni+ 1 and [wi] G 2 W , with as defined in (4.5). O 


(li ) Ijn ■ ■ ■) ) 


G 


Remark 4.7. (i) Although the set 2 W has uncountable cardinality, we have preferred this 
representation for the control actions (instead of e.g., to select W), because for each agent i 
and cell configuration 1,; the possible actions which generate a transition are bounded by the 
successor states of the agent. 

(ii) For each agent i and cell configuration lj = (i*, lj i: ..., lj N .) it holds Ur lu .] g2 wPost^(Zi; (1*, [w*])) C 
Posti(^; lj), since the transitions in TSf are associated to the specific controller selection for 
the cell configuration. < 


Next, notice that according to Definition |4.3| a well posed space-time discretization requires 
the existence of a well posed transition for each agent i. The latter reduces to the selection of an 
appropriate feedback controller for i. which also satisfies Property (P), and the requirement that 
the selected feedback controllers of the other agents also satisfy (P). Yet, it is not completely 
evident, that given an initial cell configuration and a well posed transition for each agent, it 
is possible to choose a feedback law for each agent, so that the resulting closed-loop system 
will guarantee all these well posed transitions (for all possible initial conditions in the cell 
configuration). The following proposition clarifies this point. 


Proposition 4.8. Consider system (3.1)-(3.2), let l = (l ±,..., Zjv) be an initial cell configu¬ 
ration and assume that the space-time discretization S — St is well posed, which according to 
Remark f.5 implies that for all i G Af it holds that Post^^; pr i (Q) ^ 0. Then, for every final cell 
configuration t = (l [,..., l' N ) G Posti(/i; pr 1 (Z)) x ••• x Postjv(^iv; P r jv(0)> there exist feedback 
laws 


Vi (l) {t, Xi , Xj , Xi@ , Wi ), 1 G Af , 


(4.6) 


satisfying Property (P), and w±,.. 
closed-loop system (3.1)-(3.2), (4.6 
and its i-th component satisfies 


wn £ W, such that for each i G Af, the solution of the 
(with v m = fcm,pr m (i)> wi G Af) is well defined on [0, 5t], 


Xi(5t,x( 0)) G SV,Va;(0) G D N : 
Xm{ o) = x m0 e Sl m ,m G Af. 


(4.7) 


Proof. The proof is given in the Appendix. 


□ 


The result of the following proposition guarantees that the selection of the controllers in¬ 
troduced in Definition |4.1| provides well posed solutions for the closed-loop system on the time 


interval [0,5t]. We exploit this result in Theorem 6.3 where we derive sufficient conditions for 
well posed space-time discretizations, which are also suitable for motion planning. Furthermore, 
Proposition |4.9| guarantees that the magnitude of the hybrid feedback laws does not exceed the 
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maximum allowed magnitude u ma x of the free inputs on [0, <5t], and hence, establishes consis¬ 
tency with the initial design requirement. In particular, it follows that every solution of the 


closed-loop system on [0, St] is identical to a solution of the original system (3.1)-(3.2) with the 
same initial condition and certain free input u(-) satisfying |uj(f)| < Umax, for all t > 0 and 
i £ A f. For certain technical reasons concerning the proofs in the next sections, it is convenient 
to obtain the first results of the proposition for feedback laws that only satisfy Properties (PI) 


and (P2) of Definition 4.1 


Proposition 4.9. Consider the space-time discretization S — St corresponding to the cell de¬ 
composition S of D and the time step St. Let l = (h, ■ ■. ,In) he an initial cell configuration 
and consider any feedback laws of the form 


Vi — ki^ pjv (Z) (t, Xi , Xj , XiQ , Wj), i £ Af 


(4.8) 


assigned to the agents that satisfy Properties (PI) and (P2). Then: 

(i) For each Wi £ W, i £ Af and initial condition ir(0) £ D N with Xi{ 0) = Xm £ Si i; i £ A f, the 
solution of the closed-loop system ( 3.1 )-( 3.2 (, (4.8) (with Vi = &i iPr .(j), i £ A f) is defined and 
remains in D N for all t £ [0, T), where 


T := min{(5t, min{T(xio, w^) ; * £ A/”}} 


and 


lim x(t) £ D 

t-yf- 


N 


(4.9) 


(4.10) 


Assume additionally that (P3) also holds, namely, that (P) is satisfied. Then: 

(iia) The solution x(t) of (3.1)-(3.2), ( |4.8[ ) above remains in D N foralltG [0, <$i] and satisfies 

\^i.pr i (J,)if,Xi{t),Xj(t),XiQ,Wi)\ f^U max ,Vt£ [0, St] ,i £ Af, (4.11) 


which provides the desired consistency with the design requirement (3.4) on the Uj ’s. 

(iib) There exists a piecewise continuous function v = (iq,..., vn) ■ [0, oo) -4 M. Nn satisfying 
|uj(t)| < u max , Vf > 0, i £ Af, such, that the solution x(-) above and the solution £(•) of (3.1)- 


(3.2), with the same initial condition as x(-) and input u(-), coincide on [0,<5t]. 


Proof. The proof is given in the Appendix. 


□ 


that the result of part (i) of Proposition 4.9 holds for any selection of 
, jPr .(i)(-) that satisfy Properties (PI) and (P2). Respectively, the results of 


Remark 4.10. Note 
feedback laws Uj = k. 

parts (iia) and (iib) hold for all selections of feedback laws Vi = k, pr .(i)(-) that satisfy Property 

(P)- 


In the final result of this section, we merge the results of Propositions |4.8| and T9 and show 
that each infinite discrete behaviour of the decentralized abstraction can be implemented by a 
continuous controller, which is compatible with the restrictions on the free inputs, and produces 
a continuous trajectory that satisfies the Invariance Assumption (IA). In particular, we prove 
that for each possible discrete transition sequence of the overall system which is compatible with 
the individual transition system of each agent, there exists a trajectory of the continuous time 


system (3.1)-(3.2) which satisfies (IA) and generates the discrete trajectory when sampled at 


time intervals of length St. Before proceeding, we introduce the following notion of the product 
transition system. 
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Definition 4.11. (i) Consider the space-time discretization S — St, and for each agent i £ Af, 
its individual transition system TSi as provided by Definition |4.4| The product transition 
system TSp := (Qp, Actp, — >p) is defined as follows: 

• Qp := I N (all possible cell configurations) 

• Actp := {*}[j] 

• 1 -4 P 1 ', iff l\ £ Postj(Z i; pr< (l)),V* e Af, for all 1 = (h,..., l N ), l' = (l[, ..., l' N ). 

(ii) Given an initial cell configuration 1° £ X N , a path originating from 1° in TSp , is an infinite 


sequence of states l 0 ! 1 ! 2 ... such that 1 


1 


i +1 


for all i £ N U {0}. < 


Remark 4.12. Given a well posed spac e-tim e discretizatio n S — St and an initial cell config¬ 
uration 1° £ X N , it follows from Remark 4.5 and Definition 
path l 0 ! 1 ! 2 ... in TSp originating from 1°. < 


4.11 


that there exists at least one 


We are now in position to state the result providing the consistency between the discrete 
plan level and the continuous controller implementation. 

Proposition 4.13. Assume that the space-time discretization S — St is well posed for the multi¬ 
agent system (3.1)-(3.2). Then, for each initial cell configuration f = (l®,... ,1^) and path 
fl 1 ? ... originating from f in TSp, and for each initial condition a;(0) £ D N of (3.1)-(3.2) 
satisfying Xi(0) £ S) o, i £ Af, there exists a free input v(-) satisfying |ni(f)| < v ma , x , Vf > 0, 
i £ Af, such that the solution x{t) of the system remains in D N for all t > 0 and satisfies 
Xi(mSt ) £ Sim for each m £ N and i £ Af. 


Proof. The proof is carried out by induction and is based on the results of Propositions 4.8 and 

?s of the 
Chapter 


4.9 (iib). Before stating the induction hypothesis, we provide certain basic properties of the 

Chapter 1], or 


deterministic control system (3.1)-(3.2) which can be found in 
21 - 

Recall that according to our hypotheses, the input set U of the multi-agent system consists 
of all piecewise continuous inputs v : M> 0 —> satisfying |t>j(f)| < u max , Vi > 0, i £ Af. Also, 

for each r > 0 we define the shift operator Sh r \U —>• U as 


which implies that 

v = Sh r (u) 


Sh r {v)(t) := v(t + r),\/t > 0, 


v'(t) = v(t + r), Vi > 0 


► v'(t — r)=v(t),\/t>r. (4.12) 

To the control system (3.1)-(3.2) we associate the transition map ip : A v -A D N with 

A v := {(t,t 0 ,x 0 -,v) ■■t>t o >0,x o £ D N , v £ U}, 

where ip(t, to, Xq; v) denotes the value at time i of the unique solution of (|3Tl)-([3^) with initial 
condition Xq at time io and input v(-). Notice, that by virtue of the Invariance Assumption 
(IA), (p(-) is well defined. The map ip(-) satisfies the following properties: 

• Causality. For each t > t 0 > 0, x 0 £ D N and v 1 ^ 2 £ U with v 1 |[ t0jt ) = ^ 2 |[t 0 ,t) it holds 

<p(t, to, X 0 -, v 1 ) = ip(t, t 0 , Xo\ v 2 ), 

where u 1 |[ to t ) denotes the restriction of id(-) to [to,t) (see Section 2). 


1 Following notation page 11] 
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• Semigroup Property. For each t 2 > t\ > to > 0, Xq £ D N and v £ U it holds 

V(h,t\,v{h,ta,x 0 \v)\v) = ip(t 2 ,t 0 ,x 0 -,v). 

•Time Invariance. For each r > 0, t > to > r, xq £ D N and v £U it holds 

ip(t,t 0 ,x 0 ;v) = ip(t - r,t 0 - r,x 0 ;Sh r (v)). 


Now consider an initial cell configuration 1° = (l ®,..., l^), a path l 0 ! 1 1 2 ... originating from 
1° in TS V , and an initial condition x(0) £ D N satisfying Xj(0) £ S)o, i £ Af. We will determine 
a free input v(-) satisfying \vi(t)\ < u max , Vi > 0, * £ A f, such that the corresponding solution 
x(t) := <p(t, 0, x(0);t;) remains in D N for all t > 0 and satishes Xi(mSt ) £ Iff 1 for each m £ N 
and i £ Af. The construction of v(-) is based on the following Induction Hypothesis, which 
constitutes the core of the proposition. 

Induction Hypothesis (IH). For each m £ N there exists a piecewise continuous input 
v m : R>o — > R^" satisfying |u™(t)| < u max , Vt > 0, i £ Af and such that 

v m \[ 0 , K st) = « K |[o,« ( 5 t) J VK = 1, (4.13) 

tp(t, 0, x(0); v m ) = (p(t, 0, a;(0); v K ), Vac = 1, ..., m — 1, t £ [0, nSt], (4-14) 

<Pi(mSt, 0, x(0); v m ) £ Si^yi £ M. (4-15) 


i>Proof of (IH). In order to prove (IH) for m = 1, we need to show that (4.151 is fulfi lled. 
From the fact that Sj(0) £ S)o, i £ AT and that 1° —l 1 , we deduce from Dehnition 4.11 and 
Proposition 4.8 that there exist feedback laws k i pr _( n(-) as in (4.6) (with 1 = 1°), which satisfy 


Property (Pband wi,...,wn £ W such that fl4Jl ) holds with f = 1 1 . Hence, it follows from 
iib) that there exists a piecewise continuous input v 1 : R>o —> t*" satisfying 


Proposition 
(7) | . : C rT i ax 


4.9 


, Vt > 0, i £ AT and such that 

(Pi(St, 0, s(0); v 1 ) £ 5 1 ;!, Vz £ AT 


which establishes (4.15|) for m = 1. 


Now assume that (IH) holds for certain m £ N. We will show that it is also valid for m+ 1. By 


m+1 


, we deduce from Definition 


4.11 


i,pr.(i)( - ) as in (|4.6|) (with 1 = l m ), which 


exploiting Property (4.15) of (IH) for m and that 1 
and Proposition 4.8 that there exist feedback laws k, 
satisfy Property (P), an d w ±,..., wn £ W such that ( |4.7| holds with 1' = l m+1 . Hence, it 
follows from Proposition 4.9 ’iib) that there exists a piecewise continuous input v : R>o -A R^" 
satisfying |vj(t)| < v n 


t , Mt > 0, i £ Af and such that 


<Pi(5t, 0, 0, x(0); v m )\ v) £ S^+i , Vi £ Af. 


(4.16) 


We following define v m+1 : R>o —> R 7 

v m+l {t) := 


as 


v m (t), t £ [0, mdt), 

v(t — mSt ), t £ [mdt, oo). 


(4.17) 


Then, it follows from (4.17) that v m+1 (-) satisfies (4.13) (with m := m + 1) and the latter 
implies (4.14) (with m := m+ 1) by causality. Hence, we get from (4.14) that 


ip(m6t, 0, x(0); v m+1 ) = ip(mSt , 0, x(0); v m ). 


(4.18) 
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From the semigroup property and ( |4.18 1 we deduce that 

= yj((m + l)St, mSt, tp> 

= ip((m + l)i 5t, m5t , <p(m5t, 0, x(0); u m ); v 


<p((m + l)<5f, 0, .t(0); v m+1 ) = ip((m + l)<5t, mSt , ip(m5t , 0, x(0); v m+1 ) 


v 

m+l\ 


Also, we get from (4.171 and (4.12) that 

v = Sh m5t (v m+1 ). 


(4.19) 


(4.20) 


Thus, it follows from time invariance and (4.20) that 

<p{(rn + 1 )5t, m5t, <p(m5t, 0, a:(0); u m ); v m+l ) = p(St, 0, p(mSt, 0, a"(0); u m ); Sh m a t (u m+1 )) 


= ip(5 1,0, p(mSt, 0, z(0); v m ); v). 


(4.21) 


Hence, we conclude from (4.16), (4.19) and ( 4.21| ) that ( |4.15| ) also holds (with m := m + 1) and 
the proof of (IH) is complete. <1 

In order to finish the proof of the proposition, define v : R>o —» R^" by 


v(t) := v m (t),rn £ N, t € [(m — 1 )St, m6t ), 


(4.22) 


with v m {-) as given by (IH) for each m £ N. Then, it follows from ( |4.22 1 and (IH) that 
|vi(t)| < u ma x, Vt > 0, i £ J\f and thus, by the Invariance Assumption (IA), the solution x(t) of 
the system remains in D N for all t > 0. Furthermore, it holds that 


^l[0 ,m5t) — U | [0,m<5i) 5 Vm £ N. 


(4.23) 


Indeed, for each t £ [0, mSt ) there exists k £ {1,..., m} such that t £ [(re — 1 )St, nSt). If re = to, 
then it follows from (4.22) that v(t) = v m (t). If re £ {1,... ,m — 1}, then we get from (4.22) 
that v(t) = u K (t) and thus from (4.13) that v(t) = v m (t). Hence, ( 4.23| ) is valid. Finally, from 
(4.23), (4.15) and causality we conclude that for each to £ N it holds 


Xi{m8t ) = (pi(m5t , 0, x(0); v) = ipi(m6t, 0, x(0); v m ) £ Sq 
and the proof is complete. 


□ 


We note that the result of Proposition |4.13| remains valid if we consider the product of 
the agents’ controlled transition systems as given by Definition 4.6 which will be determined 
explicitly in Section 6 for the control laws k l: \ z in (3.7). This observation is summarized in the 
following remark. 

Remark 4.14. Instead of the product TS-p formed by the agents’ individual transition systems 
TSi, consider the controlled product transition system I'.S'p := (Q' p . Aci' v .—> p ) formed 
by the controlled transition systems TS£, i £ J\f in Definition 4.6 as follows: 

• Q C V =1 N ; 

• Actj, = {*}; 

• 1 —> c v 1'. iff there exist [wi],..., [wn] £ 2 W such that l\ £ Post°(Z»; (pr^l), [iUj])), Vi € A f, for 
all 1 = (h ,..., l N ), l' = (l [,..., l' N ). 

Then, the result of Proposition 4.13 remains valid for any initial cell configuration 1° and path 
l 11 ! 1 ! 2 ... originating from 1° in TS < 
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5. Time Domain Properties of the Control Laws 


In this section we use the results of Section 4 in order to prove certain useful properties of the 
reference trajectory Xi(') an d the time domain [0,Ti(xio,Wi)) of the control laws (3.7) as spec¬ 
ified by (3.12). We proceed by providing some extra details for the dynamics as determined by 


the control law in (3.2). In particular, we assume that the ff s are globally Lipschitz functions. 
Furthermore, if we want to achieve more accurate bounds for the dynamics of the feedback 


controllers assigned to the free inputs (those will be clarified in the proof of Theorem 6.3 
the next section), we can choose (possibly) different Lipschitz constants L\,L 2 > 0 such that 


|/i(*i,Xj) - fi(xi,yj)\ <ii|(a:i,Xj) - (a^y^l, 

I fi(xi,xj) - /i(2/i,Xj)| <L 2 \{x i ,yL j ) - (yi,Xj)|, 
Vxj ,yi GD,Xj ,yj £ D Ni ,i£j\f. 


(5.1) 

(5.2) 


In order to provide some extra informal motivation on considering both constants L± and L 2 , 
we recall that in order to derive sufficient conditions for a well posed discretization, we design 
for each agent i inside a cell S). a feedback, in order to “track” a given reference trajectory (of 
i) starting in the same cell. In particular, the constant L\ provides bounds on the feedback 


term (3.9) which compensates for the deviation of agent’s i dynamics from its corresponding 


dynamics along the reference trajectory, due to the time evolution of its neighbors’ states. On 


the other hand, the constant L 2 provides bounds on the feedback term (3.11) which compensates 


for the deviation of the initial state with respect to the initial state of the reference trajectory. 

Based on the global Lipschitz assumption, we establish uniqueness of the reference trajectory 
Xi(') and provide a lower bound for the right endpoint T max of its maximal interval of existence, 
which is independent of the selection of (x it G,*.j,G) in (3.6). 


Lemma 5.1. For each tuple of reference points {xi t a, x j,g) ns in (3.6), the initial value problem 


(3.14) has a unique solution which is defined and remains in D on the right maximal interval 
). Furthermore, it holds 


[0 ,T m 


Trr 


> 


2 MLi ma,x{y/N~i : i £ J\f} 


(5.3) 


Proof. For the proof of the lemma we exploit the result of Proposition 4.9 In particular. 


we show that the solution Xi(') °f (3.14) coincides on a suitable time interval with the i-th 
component of the solution of the multi-agent system (3.1)-(3.2) under an appropriate selection 
of the initial conditions and feedback controllers for the vfs. Hence, by implicitly exploiting 
the Invariance Assumption (IA) that leads to the result of Proposition |4.9|4ia) , which is valid 


for any choice of feedback laws that satisfy Property (P), we will verify that (5.3) is fulfilled. 
In order to proceed with the proof, let x j,g) be a tuple of reference points as in 


(3.6), corresponding to a cell decomposition {S)}/ 6 z of D and a cell configuration 1, of agent 
i, and consider another cell decomposition {Sf }i 0 ex 0 of D and an initial cell configuration 
1 o 0oi> * * • > loN ) Gi To with pi'j(l 0 ) — (^oU ioji j • * •) lojN ); such that 


Ci,G 


G £? . and S? . = 


= x 


3 k ,G 


, K 1, . . . , TVj. 


(5.4) 


We have selected the auxiliary cell decomposition {Sf o }i 0 ^i 0 with the sets Sf o . consisting of a 
single element, because this slightly simplifies the subsequent analysis and also allows obtaining 
a greater (uniform) lower bound for the time T max . Next, define the time step 


5t n := 


2ML\ max{\/A : i £ A/”} 


(5.5) 
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and consider the feedback laws fc iiPr (i o ) : —>• K" given by 

k»,pr i (l o )0Cii x j) := fi(Xi,Xj,G) ~ fi(Xi,Xj) = F ith (Xi) - fi(Xi,Xj), 


with Fi \ i (•) as in (3.131 and A^ iPr< ,(i 0 ) : D Ne+1 —» K" for £ S A f\ {!} given by 

^,pr^(1 0 ) Xj(^)) • 0. 


(5.6) 


(5.7) 


Note that the feedback laws kt pr,(i 0 )( - ) for £ £ J\f\ {i} satisfy Property (P) by default. Hence, 
in order to invoke Proposition 14. 9fiia) , we show that k i pl . (i 0 )(-) also satisfies (P). Property (P3) 
is obvious, since fcj iPr .(i o )(-) is independent of t. Property (PI) follows from the corresponding 
Lipschitz property for _/)(•) and F it i.(-), since the latter satisfies the Lipschitz condition 

\Fi,ii(x) - F iM (y)\ < L 2 \x - y\,Vx,y £ D. (5.8) 


Indeed, due to (5.2) and (3.13), we have that for each x,y £ D it holds 


I Fi,u(x) - F itU (y)\ = | fi{x,x^ G ) - fi(y,x jiG )\ < L 2 \(x,x jtG ) - (y,x ilG )| = L 2 |a: - V |- 
In order to show (P2), notice that due to ( |5.5[ ) we get 

Vmax > 2M8t 0 Li\/Wi , for all i £ J\f. 


(5.9) 


Hence, we get from (4.1), (5.1), (3.5), (5.4) and ( |5.9| ) that for every Xi £ ( S° a . +H(i? max )) D D 
and Xj K £ B(xj K ,Gi Umax) nil, k = 1,..., Ni, it holds 


Ni 


\hprSo)( X i’ X j)\ =1 fi(Xi,*-j) ~ fi{Xi,XLj) | < L^Xj - X jjG | = Li ^ (Xj K - Xj^ G ) 2 

\K=1 

<L 1 ^/N~ i R max = Li^/N~iSt 0 (M + v max ) < 2M5t 0 L 1 ^/W i < u max , 


and thus (P2) holds as well, since k ipl .n g \(-) is independent of t, Xhj and Wi. Then, it follows 
from Proposition |4.9[ iia) that the solution x(t) of the closed-loop system (3.1)-(3.2), (5.6)-(5.7) 
with initial condition x(0) £ D N satisfying Xi( 0) = Xi tG , 2^(0) = Xj ljG , . .., Xj N . (0) = a gn ., G 
(and the initial state of each other agent t belonging to Sf e ) is defined and remains in D N for 
all t £ [0, St 0 }. Hence, the i-th component of the solution x(-) satisfies 


a hit) £ Dyt £ [0, <S* 0 ], 


and by virtue of ( |3.1[ )-( |3)2| and (5.6), it holds 

Xi = Fi t \ t (xi),Xi(0) = x itG ,t £ [0, J*o]. 


(5.10) 


(5.11) 


Hence, it follows from (5.11) that Xi(-) coincides with the unique solution Xi(') of (3.14) on 
[0,(5l o ] D [0,T max ), which in conjunction with (5.10) implies that Xi(t) remains in a compact 
subset of D for t £ [0, 5f 0 ] D [0, T max ). From the latter, we deduce that T max > 8t 0 . Indeed, oth¬ 
erwise Xi{t) would remain in a compact subset of D for t £ [0,T max ), contradicting maximality 
of [0,T max ). Thus, we conclude that (5.3) is satisfied. □ 


By exploiting Lemma |5.1[ it will be shown in the next section that T max is always greater 
than the maximum possible selection of the time step St for a well posed discretization. The 
latter in conjunction with the result of Lemma |5.2| below enables us to prove that in this case 
the control law ^ 3 (-) and hence also /q i ; (•) are well defined on [0, $£]. 
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Lemma 5.2. Consider a cell decomposition S of D, a time step St and select an agent i £ J\f 
and a cell configuration f = (li,lj 1 ,... ,lj N ) of i. Also, consider a tuple of reference points 
(xi'G, x j,c) os in ( |3.6| ) and the control law kij,(-) in (3.7). We assume that kij i (-) satisfies 
Properties (PI) and (P2) of Definition j.l, and that the right endpoint T max of the interval 
where the reference trajectory (3.14) is defined, satisfies T max > St. Then, for all £ Si, and 
£ W, the time Ti(xio,Wi) satisfies Ti(xio,Wi) > St, which implies that ki t i i (-) also satisfies 


Wi 


Property (P3) of Definition f.l 


Proof. Indeed, let x *o £ Si, and Wi £ W. By defining 


Xi(t) := Xi(t) 



(XiO •Ti^G),t £ [0,T max ), 


(5.12) 


with Xi(t) = Xi(t) +tWi as given in (3.15), and taking into account the definition of T(xio,Wi) 
in (3.12), we want to show that Xi (-) remains in D for more than time St. By virtue of our 
assumption that T max > St, the latter is meaningful to verify and implies that T(xio,wi) > St. 
We next show that Xi(j coincides on a suitable time interval with the i- th component of the 
solution of (3.1)-(3.2) by choosing appropriate initial conditions and feedback laws that satisfy 


(PI) and (P2). 

Let Xio £ Si,, Wi £ W, consider an arbitrary initial cell configuration 1 with pr,(l) = 1,;, 
1 = (l i,..., In), and assign the feedback law k i pr .^ = kn, (as the latter is given by (3.7)) to 
i and the feedback laws A^ lPr (i) := 0 to the rest of the agents i £ Af\ {*}. It also follows from 
the assumptions of the lemma for i, and trivially for the other agents, that the feedback laws 
satisfy Properties (PI) and (P2). Thus, we can use the result of Proposition |4~9{ i). By selecting 
an initial condition x(0) £ D N with x,(0) = x,;o and Xj m (0) £ Si m , m = 1,..., iV„ and recalling 
that Wi £ W, we get from Proposition |4.9[i) that the 7-th component of the solution satisfies 


Xi(t) £ D,\/t £ [0, T), T := min {St,T(x i0 ,Wi)} 

lim Xi(t) £ D. 
t-*f- 


(5.13) 

(5.14) 


We proceed by showing that Xi(t) 
equivalently, that 


Xiif), for all t £ [0,T), with T as given in (5.13), or 


Xi(t) 


Xi(t) + tWi + 



{x i0 - x ijG ),Vt £ [0, T) 


(5.15) 


Indeed, from (3.15), (3.14), (3.1)-(3.2), (3.7), (3.9) and (3.13) we have that 


ii(t) = F i>u (xi(t)) + Wi, 

Xi(t) = F iiU (xi(t)) + k it i it2 (xio) + ki'i i>3 (t-,Xio,Wi). 


By recalling that x,(0) = ar»(0) = x^ and that due to (3.12) and (5.13) it holds T < 

T(xio,Wi) < T max , and thus Xi(')> Xj(-) and u, 3 (-) are well defined on [0,T), it follows from 
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(3.10), (3.11) and (3.15) that 


Xi(t) - Xi(t) = x i0 - x it G + / [Fi,u{xi(s)) - 

Jo 

+ h,li, 2 (Xio) -I- XiO,Wi) - Wi]ds 

{ x i0 ~ x i,c) + J [-Fi,!<(£*(«)) 


= 1 - 


- F, 


,L (xi(s) + (i - (x l0 - x iiG )) ds,\/t G [0, T). 

Hence, we get from (|5.8[) that for all t G [0,T) it holds |xj(i) — Xj(t) — (l — ^) (xio — x^g)I < 


f 0 L 2 |x,;(s) — x*(s) — (l — fi) (xjo — x^g)) ds. Application of the Gronwall Lemma, (3.15), and 


the fact that T < T max , imply that (5.15) holds. 


We are now in position to prove that T(x,;o, w,) > St. Indeed, suppose on the contrary that 
T(xio,Wi) < St , which by virtue of the assumption that T max > St, and (5.13), implies that 
T(xio,Wi) < T m ax and T = T(xio,Wi). From the latter, together with ( |5.12 1, (5.13), (5.14) and 
continuity of Xj(- ), we get that x l (T(x i0 ,»,;)) = ®»(*) = x, ; (t) G 

D. Hence, from (5.12), the deduction that T(xio,u>i) < T max and continuity of Xi(-), it follows 
that there exists £ > 0 such that Xj(t) G D for t G [T(xio, Wi), T(xi 0 , Wi)+e), which contradicts 
(3.12). Thus we conclude that T(xio,Wi) > St, which establishes validity of (P3). □ 


6. Well Posed Space-Time Discretizations with Motion Planning Capabilities 


In this section, we exploit the controllers introduced in (3.7) to provide sufficient conditions 


for well posed space-time discretizations. By exploiting the result of Proposition |T8] this frame¬ 
work can be applied for motion planning, by specifying different possibilities for transitions 
for each agent through modifying its controller. Consider again the system (|3T|)-(|T2|) , a cell 
decomposition S = {Si}i^x of D and a time step St. In addition, consider the least upper 
bound on the diameter of the cells in S , namely, 


dmax := sup{diam(5'i), l G I}, 


( 6 . 1 ) 


which due to Definition [3T] is well defined. We will call d max the diameter of the cell decom¬ 
position. Our goal is to determine sufficient conditions relating the Lipschitz constants L 1; L 2 ; 
the bounds M, u max for the system’s dynamics, as well as the space and time scales d max and 
St of the space-time discretization S — St, which guarantee that S — St is well posed. According 


to Definition 4.3 establishment of a well posed discretization is based on the selection of ap¬ 
propriate feedback laws which guarantee well posed transitions for all agents and their possible 
cell configurations. For each agent i G AT and cell configuration 1, ; = {h,lj 1 ,... ,lj N .) of i let 


(xj,G) x j,g) be a reference point as in (3.6). We consider the family of feedback laws given 


in ([T9]), pTOl), (|3.11[), and parameterized by x^ G Si, and Wi G W. The function Fi.pO) is 


given in (3.13), and Xi(') is the reference solution of the initial value problem (3.14), defined 
on [0, T max ). Recall that the parameter A in (3.8) provides the portion of the free input that is 


exploited for planning. Thus, it can be regarded as a measure for the degree of control freedom 
that is chosen for the abstraction. In the following results, we also introduce an additional pa¬ 
rameter /i which provides a lower bound on the minimum number (> 1) of discrete transitions 
that are possible from each initial cell, as will be clarified in the corollary at the end of the 
section. Before proceeding to the desired sufficient conditions for well posed discretizations and 
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their reachability properties, we prove the auxiliary Propositions |6.1| and |6.2| Proposition |6.1| 
below provides bounds on the hybrid control laws fc* i 4 (-) in (3.71. 

Proposition 6.1. Consider a cell decomposition S of D with diameter d max and a time step 
St. Also, for each agent i £ Af and cell configuration k = (h, / Ji; • • ■, lj N .) of i let (x^g, x j,g) be 
a reference point as in (3.6 1 and consider the feedback law (•) in (3.71. Then its components 
ki,k, i(')> ki,k, 2(0 an d ki,k, 3(') as given in (3.9|, (3.101 and (3.11), respectively, satisfy the 
bounds 


| ki,k ,1 x j ) I — T din ax ), V.X‘ ? ; £ D, 

x jm £ ( Si m + B(R max )) n D,m = 1 ,... ,N t , 

\ki,k, 2 (xio)\ — ^d max ,Voqo £ S[ i , 

XiO, wf) | ^ ((5tAu max -p rf max ) T Au max , 

vt G [0, <5i] n [0, T(x i0 , Wi )), a; i0 € S Zi , to* G W. 


( 6 . 2 ) 

(6.3) 

(6.4) 


with i? max as given m (4.1). 


Proof. Indeed, in order to show (6.2) let Xj G D Ni satisfying Xj m £ (S Zm + B(R max )) C\D,m = 
1 ,..., IVj. Then, for each m = 1,..., IV* there exists Xj m with Xj m £ iS) m and — Xj m \ < R max . 
Hence, from the latter together with (3.9) and (5.1), we get 

|fci,l 4 ,i(a;i,Xj)| < L 1 \{x h - x juG , ■ ■ -,x jN - Xj G )I 


< L\ 


( Ni 

E( 

\m=l 
f Ni 


< U E (■^max H - ^max ) J — k M, (/fmax -p d max ), 


which establishes (6.2). Furthermore, by recalling that it follows directly from (3.10) 

that |fcj ) i ii 2 (a:io)| = — x i,c\ and hence, that (6.31 is satisfied. Finally, for ki } l^O) we get 

from (3.11) and ( |5.8| ) that 

IhxM^Xio, «j*)| < L 2 |(xi(t) + tWi + (l - ji) (xio - x itG )) - Xi(f) | + K|, 


which due to (3.8) implies validity of (|6.4|). 


□ 


Based on the result of Proposition |6.1| we next provide conditions on d max and St which 
guarantee that the feedback laws &i,ii(') satisfy Property (P). Additionally it is shown that the 


radius r introduced in (3.16) satisfies a design requirement which is related later in Corollary 6.7 
to a lower bound on the number of possible transitions through the parameter p. 

Proposition 6.2. Consider a cell decomposition S of D with diameter d max , a time step St, 
the parameters A £ (0,1), p > 0 and define 


L := max{3L 2 + 4Li i £ J\f}, 


(6.5) 


with Li and L 2 as given in (5.1) and ( |5.2[ ). We assume that A, p, d max and St satisfy the 
following restrictions, as provided by the three cases below: 
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Case I. 0 < p < 

^max ^ ( 0, 


(1-A) 2 <, 


4 ML 


( 6 . 6 ) 


St £ 


(1 - A)tJ max - a/(1 - A) 2 i4 ax - 4MLd max (1 - A)v max + ^/(l - A) 2 z^ ax - 4 MLd„ 


2 ML 


2 ML 


Case II. AW < p < 


4A 


^rnax £ 6, 


2(A(1 — \)fi — 2A 2 )r 


<5* e 


2Xv n 


li 2 ML 

(1 - A)v max + \J(l- A) 2 i>^ ax - 4 M~Ld~ n 
2 ML 


(6.7) ' 

( 6 . 8 ) 

(6.9) 


2(A(l-A)M-2A 2 )< ax (1 — A) 2 


H 2 ML 


AML 


( 6 . 10 ) 


and St satisfies (6.7). 


4A 


and St satisfy (6.8) and (6.9), respectively. 


Case III. p > 

Then, the intervals in Cases I, II, III are well defined, and for each agent i £ A/", cell configu¬ 
ration li = (li, lj 1 ,..., lj N ) of i and reference point (xi : c, x j,c) as in (3.6) the solution Xi(t) °f 


(3.14) is defined and remains in D for all t £ [0,5t]. In addition the feedback law ki^(-) in (3.7) 
satisfies property (P) and the distance r as defined in (3.16) satisfies the design requirement 


r > —d 

1 — 2 Uma x 


( 6 . 11 ) 


Proof. The proof of the fact that the intervals in Cases I, II, III are well defined is provided in 

that for any reference point (xi t G,Xj,G) as in 


the Appendix. Also, it follows from Lemma 


5.1 


(3.6) the solution Xi(t) of (3.14) is defined and remains in D for all t. £ [0,2),; 
of (5.3) and the assumed bounds on St in Cases I, II, III, that 


> St, 


*), and by virtue 


( 6 . 12 ) 


[0, <5£]. We break the subsequent proof in the 


establishing thus that %j(t) £ D for all t £ 
following steps. 

STEP 1: Verification of Properties (PI) and (P2) for the feedback law (3.7) for 
dm ax — St as given by Cases I, II, III, in conjunction with the design requirement 


(6.11). In this step we prove that the proposed feedback law ( |3J| satisfies Properties (PI) and 
(P2). Verification of (PI) is straightforward. Thus, we proceed to show that (4.2) holds, which 
implies (P2), and simultaneously, that ( 6.1 1| ) is fulfilled. By taking into account (3.7) and the 
result of Proposition 16.ll namely, (16.21), (16.3) and (16.41), we need to prove that 


L iy /Ni(R max + d m; 

+L 2 (StXv n 


x) T ^ i d max 

ix T d max ) -f- A e Tn ,, x £ v ,, 


(6.13) 


By recalling (4.1), (3.5) and the fact that A £ (0,1) we get that StXv max < . Also, from the 

fact that d max and St are selected according to the Cases I, II, III, it follows from elementary 
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calculations which are provided in the Appendix that 

dmax 111 Am a x- 


Hence, it suffices instead of (6.131 to show that (2 Li^/N~i + | L 2 )R„ 
which by virtue of (4.1) is equivalent to 


id 

St “max 


(M + v max )(2Li\/Ni + -L 2 )St 2 — (1 — A )v max 5t + d max < 0. 


(6.14) 
i (1 A)ll max , 

(6.15) 


By taking into account (3.5), it suffices instead of ( 6.15| ) to show that M(3L2 + ALi^/TSTfSt 2 — 
(1 — X)v max St + d max < 0 which by virtue of (6.5) follows from 


ML5t 2 - (1 - A )v max St + d max < 0. 

In order for the above equation to have real roots, it is required that 


(1 — A) 


V. 


- 4 MLd max > 0 


(1-A) 2 i4 

AML 


(6.16) 


(6.17) 


Hence, by collecting our requirements (6.17), (6.16), (6.14) and (6.11) together with the fact 
that d max > 0 we have 


0 <C d max < 
(1 A)z; maj 


(l^A) 2 <ax 
4 ML 

- V / (l-A) 2 < ax -4 MLd, 


2 ML 


max , e , (1 A)l? max 
< Ot < - 


(6.18) 

+ Vi 1 ~ A) 2 u^ ax - 4MLd max 


2 ML 


M- 


2Xv n 


< 5t, 


< St. 


(6.19) 

( 6 . 20 ) 
( 6 . 21 ) 


We can then show that for all Cases I, II and III as in the statement of the proposition the 
above requirements are satisfied and hence, that (P2) holds. The proof of this fact can be found 
in the Appendix. 

STEP 2: Verification of Property (P3). In order to show (P3), it suffices to prove that 
for the given selection of A G (0,1), fi > 0, d max and St as provided by Cases I, II, III, the agent 
i and the cell configuration 1; it holds T(xio,Wi) > St, for all Xw £ and Wi £ W. The latter 
is a direct consequence of (|6.12 ) and Lemma 5.2 


□ 


We are now in position to state our main result on sufficient conditions for well posed 
abstractions. 

Theorem 6.3. Consider a cell decomposition S of D with diameter d max , a time step St, the 
parameters A £ (0,1), p, > 0 and assume that X, pt, d max and St satisfy the restrictions of 
Proposition 6.2. Then , for each agent i £ AT and cell configuration k = (fi, lj 1 ,..., lj N ) of i 
we have Post7[q; If) 0, namely, the space-time discretization is well posed for the multi-agent 
system (3.1)-(3.2). In particular, for any tuple of reference points (xi t c, x j,c) as in (3.6) and 
corresponding reference trajectory Xi(-) of i as given by (3.14) it holds 

B{Xi(St)',r) C D, (6.22) 

Post iik; k) 3 {l £ T : S t C B{ Xi {8f)-, r) ^ 0}, (6.23) 


where r is defined in (3.16). 
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Proof. For the proof, pick i £ A f, 1 , = (k, lj lt . .., lj N .), (x^g, x j,g) as in (|-i.6[) and notice that 


by virtue of Proposition 6.2 the reference trajectory Xi(') is well defined on [0, St]. In addition, 
consider the control law in (3.71. Then, it follows again from Proposition 

latter satisfies Property (P). Next, for each x £ B(xi{St);r) define 


6.2 


that the 


Wi(= Wi(x)) := 


x - xSt) 

St 


Then, we get from (3.16) that «y| < T = Au max and hence, by virtue of (3.8) that 

Wi(= Wi(x)) £Wyx £ B(xi(St)\r). 


(6.24) 


(6.25) 


In order to prove the theorem, we need to verify that (6.22) and (6.23) are fulfilled. 


Proof of (6.22). In order to show (6.221, pick x £ B(xi(St);r ), Wi as in (6.24) and recall 


that the control law fci.pO) satisfies Property (P). Then we have from (6.251 that Wi £ W and 
thus, we get from Property (P3) applied with aq 0 = Xi g and the selected parameter Wj that 
T(xi’G,u)i) > St. From the latter and (3.12) we obtain that X;(<5t) +WiSt £ D, which by virtue 


of (6.24) implies that x £ D and establishes validity of (6.22). 


Proof of (6.23). For the verification of (6.23) it suffices to prove the following claim. 


Claim II. Consider the control law ki t i 4 (-) above and pick any Wi £ W. Then, for any initial 
cell configuration 1 with pr 7 -(1) = l i; 1 = (l lt ... ,In), £ £ Af\{i} and selection of feedback laws 
in (4.4) which satisfy (P) the following hold. The solution of the closed-loop system (3.1)-(3.2 ), 
(3.7), (4.3), (4.4), with the selected parameter Wi for fcj,!*(•), is well defined on [0,(5f] and satisfies 


Xi(5t,)(:= Xi(8t,x( 0))) = Xi(St) + Stw t = Xi(St), 


(6.26) 


for all x(O) £ D N with x m o £ Si m , m £ Af and w m £ W, m £ Af\ {*}, with the last equality 
in (6.26) being a consequence of ( 3.15| ) (see also Fig. 2 in Section 3). < 

Indeed, let any l £ X such that 


Si C B(xi(St); r) £ 0. 


(6.27) 


In order to show that l £ Post i(li;h), i.e., that the transition /, —->■ l is well posed, it suffices 
according to Definition |4.3| (i) to verify that there exists Wi £ W such that Condition (C) holds 
with the control law fci,p(-) above. By exploiting (6.27), we pick x £ Si fl B(xi(St); r) and Wi in 
(6.24) as the parameter for ki t i 4 (-), which by virtue of (6.25) satisfies wi £ W. Thus, it follows 
that the conclusion of Claim II is fulfilled with (6.24| and (6.26) implying that Xi(St) = x £ Si. 


Hence, Condition (C) is satisfied and we conclude that li -4- l is well posed. It thus remains 
to verify Claim II. 

Proof of Claim II. Let x, fcj ^ (•) and Wi(= Wi(x)) as in the statement of Claim II. We 
first note that due to Proposition |4.9[ iia), the solution of the closed-loop system is defined 
and remains in D N on the whole interval [0,<5t]. In order to show that Xi(St) = Xi(St), we 
show that Xi(-) is an appropriate modification of the trajectory Xj(-). In particular, it holds 
Xi(t) = Xi(t ) + (l — Jj) (xio — Xi, G ),Vt £ [0,<5f], which implies the desired result. The proof of 
the latter is based precisely on the same arguments used for the proof of (5.15) in Lemma 5.2 


and is therefore omitted. Hence, we conclude that Xi(St) = Xi(5t) and the proof is complete. □ 


Remark 6.4. Notice that the reachable cells provided by the left hand side of (6.23) are a subset 


of the reachable cells from the specific cell configuration of i. In addition, these cells depend on 
the reference points (xi : c, x-j g), since a different selection will correspond to a different control 
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law ki^ t and will lead in principle to another reference trajectory, thus, modifying the center of 
the ball B(xi(dt);r). <s 

Based on the above remark, we next show that by selecting for each agent i and cell config¬ 


uration I, a tuple of reference points (cc^Gi x j,g) and the control law (3.7), the right hand side 
of (6.23) provides the set Ur w .] e2 wPostf(ij; (1*, [wjj])) of the controlled transition systems TSf 
in Definition [TGI 


Corollary 6.5. Consider a space-time discretization S—St satisfying the hypotheses of Theorem 
\6.S\ Also, select for each agent i and cell configuration k, a tuple of reference points (x^g , x j,g) 
as in (3.6), the control law ki^{-) in (3.7) and consider the corresponding reference trajectory 
Xi{') given by (3.14). Then, for each i, f = (l i ,lj 1 ,... ,lj N f) and l £ T the actions in 


(4.5) for the specification of the controlled transition system TSf are given as 

x~Xi(8t) 


M(w) = 


5t 


: X £ Si D B{xi{6t)\r) 


with r as in (3.16), and it holds 

u K]e2^ p ost-(Zi; (Z*, [w;])) 

= {lel:S l nB( Xi (6ty,r)^(/)}. 


(6.28) 


(6.29) 


Proof. Proof of (6.28). We first consider the case where Si fl B{xi{8t)]r) ^ 0 and show that 
[wi](L ,n D { X ~ X st '■ x e SinB (xi(6t);r)}. Indeed, let x £ Si r\B(xi(St)-, r) and Wi as given by 
(6.24), which by virtue of (6.25) satisfies w, £ W. Thus, from Claim II in the proof of Theorem 
we obtain that I,;, fci,i;(‘)> w i and l satisfy Condition (C), since (6.24) and (6.26) imply that 


6.3 


Xi(St) = x £ Si. Thus, we obtain from (4.5) that - (i is i)• In order to prove the 

reverse inclusion, let w, £ C W and assume on the contrary that 


x - Xi($t) 

St 


■ x £ Si n B(xi(St)',r) 


(6.30) 


Let x = Xi(bt) + StWi and notice that due to ( |3.8[ ) and (3.16) it holds x £ B(xi{St);r). In 
addition, by exploiting Claim II applied with this selection of w^, we obtain from (6.26) that 
Xi(St) = Xi(fit) + dtwi = x. Also, since Wi £ ['U'ijp,,;) we get from (4.51, namely, the fact that 
1 1 , fci,h(')) w h l satisfy Condition (C), that Xi(St) £ Si. Thus, x £ Si fl B(xi(St)-,r) which 
contradicts (6.30|). 


Finally, in order to verify (6.28) for the general case we need to show that [iDi](i^,z> = 0 when 
Si fl B(xi(St)\r) = 0. Notice that the latter is equivalently written as 


\x - Xi($t )| > r,\/x £ Si. 


(6.31) 


Hence, suppose on the contrary that (6.31) holds and that [u , i](i il z) ^ 0. Then, by picking any 
Wi £ ['Wj](i ii q C W , we obtain from (6.26) in Claim II that Xi(St) = Xi(^t) + dtwi. On the 
other hand, since Wi £ [wi](L,i), we get from (4.5) that Xi(St) £ Si. Hence, it follows that 
Xi($t) + StWi £ Si which by virtue of (6.31) implies that |ui,| > A. Thus, we obtain from (3.8) 


and (3.16) that Wi ^ W, which is a contradiction. 


Proof of (6.29). This follows directly from (6.28) and the definition of transitions in TSf as 

' □ 


provided by Definition |4.6| 


Remark 6.6. It is noted that if we select W in (3.8) as the open ball int(B(An max )) we will 
obtain in (6.23) the cells which have nonempty intersection with the open ball int(7?(xi(<5t); r)). 
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In this case, it follows from the properties of a cell decomposition that each transition system 
TSf is deterministic. The latter relies on the requirement that the interiors of the cells in the 
decomposition are disjoint and the fact that when the intersection of a cell with this ball is 
nonempty, it will also have nonempty interior. < 


The results of Theorem |6.3| and Corollary |6.5| can be utilized for motion planning tasks 
through the following procedure: 

Step 1. Given the Lipschitz constants L\, L 2 and the bounds M, v max on the agents’ dynamics, 
pick design parameters A, g and select a well posed space-time discretization S — St for the 
multi-agent system based on Theorem |6.3| 

Step 2. Fix a reference point for each cell Si of the decomposition {S)}; 6 x- Then, derive 
the controlled transition system TSf of each agent i as follows. For each cell configuration 
1,; = ... ,lj N ) compute the endpoint x,(c>f) of the reference trajectory (3.14) at time St, 

corresponding to the reference points Xi } c, Xj, ,c:-. • • ■, Xj N . ,g of the cells Si z , ,..., Si.. ^ , as 

selected at the beginning of Step 2. Then, specify the cells which have nonempty intersection 

with B(xi{St);r), in order to obtain all the transitions k — l[ to the cells in (6.29), with 


[iUj] = as given by (6.28). Also, note that the actions [w^] do not need to be specified 

until Step 4. 

Step 3. Find a path l 0 ! 1 ! 2 • • • in the controlled product transition system TS£, defined in 
Remark 14.141 which satisfies the plan and project it for each agent i to a sequence of transitions 
(iJ.M 1 ) 

/0 _vc /l _v c /2 . . . 

Step 4. Select the control laws to implement the individual transitions by the continuous time 

(ir.Kn 

system as follows. For each transition l™ —pick any parameter Wi € [ Wi] m = 
; m+ 1 ,, with the latter as given in (6.28), and apply the control law (3.7) with the selected 
Wi- 

The following corollary provides a lower bound for the minimum number of cells each agent 
can reach in time St, depending on the selection of the design parameter ^ for the space-time 
discretization. 

Corollary 6.7. Consider a cell decomposition S of D with diameter d ma x, a time step St, and 


parameters A £ (0,1) , p > 0 such that the hypotheses of Theorem 6.3 are fulfilled. Then for 
each agent i £ N and each cell configuration of i, there exist at least \_p n \ + 1, if gT N, or 
L/z n J, if g n £ N possible discrete transitions. 


Proof. In order to prove the result, we need by virtue of (|6.23|) to show that 


lT n \ 

[ h n l 


1, 


if g n $ N, 
if p n £ N, 


(6.32) 


where ff denotes the cardinality of a set. By using the notation Vol(S') for the volume (Lebesgue 
measure) of a measurable set S C 
that for each l £l it holds 

dmax A A f dmax 

~ 2 / J = 


Vol (Si) < Vol B 


it follows from (6.1) and the iso-diametric inequality 
0(n) '■= S max , (6.33) 


where 


P(n) :=Vol(B(l)), 
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namely, the volume of the ball with center 0 and radius 1 in M ra . It then follows from (6.22) 


namely, that B(xi(St);r) C D, (6.33) and the fact that due to Definition 3.1 it holds U i^xSi = 
D , that 

{ I Vol {B(Xi{St)-,r)) I , ,r ( B(xi(St);r )) i j^r 

#{l £l: Si ± 0} > j j ’ , f (B(£%$t)-,r)) c ^ ( 6 - 34 ) 


By taking into account (6.11|) and (6.33), we get that 


Vol (B( X i(Sty,r)) ^ (fd m ax) n /3(n) _ 

^ m - fJj 


s„ 


— Id 


/3(n) 


(6.35) 


and thus, (6.32) is a direct consequence of (6.34[) and (6.35). The proof is now complete. □ 


Finally, we provide certain upper bounds on the complexity of the controlled transition 
system of each agent as a function of A, namely, the part of the input that is exploited for 
reachability purposes. Therefore, given a finite cell decomposition S = {Sz}; e x of a bounded 
domain D , it is convenient to introduce the length 


d in = 2sup{I? > 0 :\/l £ l,3x £ Si,B(x , R) C Si} 


(6.36) 


corresponding to the maximum diameter of a ball that can be inscribed in all cells. The following 
corollary provides the corresponding complexity result. 

Corollary 6.8. Consider a bounded domain D admitting a finite cell decomposition S of diam¬ 
eter d max , a time step St and a parameter A £ (0,1), such that the hypotheses of Theorem 6.3 
are fulfilled with p = 0, implying that d max and St satisfy (6.6) and (6.7), respectively. If in 
addition d max 


is the maximum possible diameter that satisfies (6.6), i.e., d max = *— .L**' , 
then the cardinalities of the state set Q\ and transition relation —of agent ’s i individual con¬ 
trolled transition system are upper bounded by C\ and C 2 ^ 1 _ A ^ 2 (K i +i)+i)r. > respectively, 

with Ci = (cv l'f ) > C 2 = C^ i+1 ( 4 )”, c = d in as given in ( |6.36| ) and Vol (D), 


Vol(I3( 4 )) being the volume of the domain D and the ball with radius 4 in R", respectively. Fi¬ 


nally, the cardinality of the state set Q%, and transition relation 


of the controlled product 


transition system are upper bounded by C[ ^ 1 _ A L pj V „ and C 2 ( 1 _ a L ) 3 jv^ , respectively, with C[ = Cf 
and C 2 = Ci 


N 


Proof. Notice first, that when d max = ^ V 


4 ML 


St 


t follows from (6.7) that necessarily 

(1 A) U max 

2 ML ‘ 


(6.37) 


Next, based on the finiteness hypothesis of S which implies that the length dj n as defined in 


(6.36) is positive, we provide an upper bound on the number of cells in S and the number of 


cells that may intersect any reachable ball in D with radius r, as the latter is given by (3.16). 


First, notice that due to (6.36) it holds 


Vol(Sj) < Vol(B(^)),VZ G J. 


(6.38) 


In addition, in order to obtain a bound on the cardinality of X, we exploit the facts that the cells 
of the decomposition cover D and have disjoint interiors, which implies that for any subset of 
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cells, the volume of their union equals the sum of their individual volumes. The latter properties 
in conjunction with (6.38) imply that 

Vol(Sj) = Vol(D) =► Vol(B(^)) < Vol(D) =► 


Vol(U ie i5 ; ) = Vol(D) = 
Vol(-D) 


lex 


m< 


Vol (D) 1 

Vol(S(^)) “ Vol(B(i))< 
Vol(D) f AML 


lex 
Vol {D) 


Vol(S(|)) (cd max ) r 


Vol(B(|)) \c(1 — A) 


Vol(D) / 4ML \ n 
Vol(B(i)) 0"-A) 


2 n 


= c\ 


(l-A) 


In ' 


(6.39) 


with c and Ci as given in the statement of the corollary. We next proceed to determine an 
upper bound on the number of cells which intersect any reachable ball B(xi{dt)',r). Note first, 
that by the definition of d max it follows that 

{l G 1: Si n B(xi(6t); r) £ 0} C {l G I: Si C B(xi{dt); r + d max ) ^ 0} => 

#{/ € 1: 5; n B(xi{St); r) 0} < #{Z € Z : Si C B{xi{St)\r + d max )}. (6.40) 


In addition, by taking into account as above that the cells of the decomposition have disjoint 
interiors, we get that 

Vol ([J{Si G S : Si C B(xi{dt);r + d ma x)) C Vol(.B(xi(<ft); r + d max )) => 

^ Vol(S)) < Vol (B(xi(6t); r + d max )) => 

{ieI:SiCB(xi(<5t);i-+d max )} 

Vol(B(xi(<5t); r + d max )) 


#{Z £l: Si C B(xi(St);r + d max )} < 


Vol(£(^)) 


[H, 

dj n 

, 2 


(6.41) 


From the selection of d max in the statement of the corollary, St in (6.37) and the definition of 
r we get that 


r + d B 


Xv 


(1 —A)u m; 
max 2ML 


d-A) 2 ^„ x 

AML 


(l-A) 


Cdrr 


< 


2ML maX 

ca max 


ML 

AML 


(l-A )- 


(6.42) 


Due to (6.39) it follows directly that the cardinality of the state set of each agent’s individ¬ 


ual (controlled) transition system is upper bounded by C\ ■ I n order to obtain the 

corresponding bound for the transition relation, note that due to (6.28) and (6.29), for each 
1* = {hi Iji , ■ ■ -, lj N .) it holds that 


k 


(ii>N) 


l'i iff l'i € {l G Z : Si n B(xi(St)-,r) ^ 0} and [ry<] = K>] (!<,*') 


(6.43) 


and that each action is uniquely determined by the successor cell l\. Thus, the cardi¬ 

nality of each agent’s i transition relation is evaluated by summing the numbers of the possible 
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of (6.43) that 


successor cells over all the possible cell configurations of i. This observation implies by virtue 

E #{ieJ:S«nB(xi(<S*);r) ? 40}. (6.44) 

Lex 2 ^ 1 


Hence, by exploiting ( |6.39| ), ( |6.40[ ), ( |6.41[ ) and ( |6.42[ ), we get that 


# 


< 


< c 


(1-A) 


2n 


Ni+1 


= C{ 


Ni+1 


(l-xy 


= Co 


-a y 


i 


(6.45) 


c J (1 - A)( 2 ( Ari+1 )+ 1 ) 71 ^ 2 (1 - A)( 2 ( iV >+ 1 )+i)^ 

with C 2 as given in the statement of the corollary. Finally, from the definition of the product 
controlled transition system, we obtain from (6.391 that the cardinality of its state set satisfies 

1 


#Qv = = mr < o 1 


(1 ~ A) 2 


N 


= ci 


N 


1 


(l-A) 


2 Nn 


= c[ 


1 


(l-A) 


2 Nn ' 


with C[ as given in the statement of the corollary. In addition, we obtain that the cardinality 
of the transition relation is given by 

N 

lex™ i=i 

where for each 1 G 1 N and i G Af, Xi,pr (l)(St) denotes the reference trajectory corresponding to 
the cell configuration pr ; -(1) of agent i. Thus, we get from (6.39), (6.40), (6.41) and (6.42) that 


N 


# 


<#x"n - 


i= 1 


(l-A)' 


<m) 


N 


Nn 


1 


< C\ 


(l-A) 2 


N /^\ Nn 


(1 - A)^" 

Nn 


(l-A) 


Nn 


= cf 


N 


1 _-c* 1 

(1 — X) 3Nn 2 (1 — X) 3Nn ’ 


with C' 2 as given in the statement of the corollary. The proof is now complete. 


□ 


7 . Example and Simulation Results 

As an illustrative example we consider a system of four agents with states X \, £2, £3, £4 G R 2 , 
whose initial conditions lie inside the circular domain int(R(i?))(= {x G R 2 : |x| < R}) with 
center zero and radius R > 0 . Their dynamics are given as: 

£1 = sat p (x 2 - £1) + s(xi) + i>i, 

£2 = g{x 2) + v 2 , 

£3 = sat p (£ 2 - £3) + g(x 3 ) + v 3 , 

£4 = sat p (£ 3 — £4) + p(£ 4 ) + u 4 , ( 7 . 1 ) 

where the function sat p : R 2 —> R 2 is defined as sat p (x) := x if |x| < p; satp(x) := j^£, if 
\x\ > p. The agents’ neighbors’ sets in this example are A/i = { 2 }, J\f 2 = 0, A /3 = { 2 }, A /4 = {3} 
and specify the corresponding network topology. The constant p > 0 in satisfies p < R 
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and represents a bound on the distance between agents 1, 2, and agents 2, 3, that we will require 
the system to satisfy during its evolution. The function g(-) is defined as 

if |rc| < R — 

M)r> if R-$<\x\<R, (7.2) 

if R < |*| 

and determines for each agent a repulsive vector field from the boundary of int(B(i?)) when 
the agent is located in intin order to ensure invariance of the agents’ trajectories inside 
this circular domain. 

We next show, that if the initial distances between agents 1 and 2 (and similarly for agents 2 and 
3) is less than p, it will also remain less than p for all positive times, for an appropriate bound 
on the magnitude of the free input terms V{. By selecting the energy function H(xi,x 2 ) '■= 
\ |*i — X 2 1 2 and evaluating its derivative along the right hand side of ( |7.1[ ), we obtain that 

V = (xi - x 2 , x 2 X\ T g(x i) + vi - g(x 2 ) - v 2 ) 

< -|*i - x 2 \ 2 + (*i - x 2 ,g(x i) - g(x 2 )) + 2|xi - x 2 |u ma x,if \%i - x 2 \ < p (7.3) 

V = {x\ - x 2 , t- - -r(x 2 - xi) + g(x i) + V! - g(x 2 ) - v 2 ) 

\x 2 - Xi\ 

< ~(p - 2v max )|*i - * 2 1 + (Xi - x 2 ,g(xi) - g(x 2 )), if \xi - x 2 \ > p (7.4) 


r 

g(x) := < {(R - f) - 

P X 

\ 2 |*| • 


where (•,•) denotes the inner product in R 2 . Next, notice that 

(x - y, g{x) - g{y)) < 0, Vx, y € R 2 


(7.5) 


Indeed, assume that without any loss of generality it holds |x| > \y\. Then, it follows from 
(7.2) that there exist a > /3 > 0, such that g{x) = —ax and g(y) = —fly. Hence, we get 


that (x - y,g{x) - g{y)) = (x - y,-ax - /3y) < -a|x | 2 + (a + /3)\x\\y\ - /3\y\ 2 . By evaluating 
the discriminant of the latter second order expression, we obtain that it is always nonpositive, 
which implies (|7.5|). By additionally assuming that 


‘'max 


P 


(7.6) 


we obtain from (7.3), 
|xi( 0 ) - x 2 


7.4[ ) and (7.5) that V < 0 when |xi — x 2 | > p. Thus, it follows that 
< p implies that |xi(£) — x 2 (t)\ < p for all positive times. Analogously, by 


considering the function H(x 2 ,x 3 ) := ||x 2 — x 3 \ 2 it follows that the same holds for | x 2 (t) — 
x 3 (t)|. Furthermore, under the selection of x max = |, it can be deduced (along the lines of the 
corresponding result in El) that the circular domain remains invariant for the dynamics of the 
system. Finally, we obtain from (7.1) and (7.2) the following dynamics bounds and Lipschitz 
constants in (3.3), (5.1) and (5.2), respectively: 



(7.7) 


Thus, it follows that system (7.1) satisfies all requirements for the derivation of well posed 
discretizations. 


In this example, it is also assumed that the reference point of each cell of the square partition 
is the center of the square. This enables us to obtain the following improved bounds on the 


















DECENTRALIZED ABSTRACTIONS FOR MULTI-AGENT SYSTEMS UNDER COUPLED CONSTRAINTS29 


feedback laws in (3.7), for their corresponding values of t, x. t , Xj and wp. 


T ^-j ) I — Li | v max ^)St -\- 


2 St 


\ki,U,3(t\Xio,w)\ < L 2 ( Xv 

max St + 


Xv n 


Thus, in order to verify Property (P2), we need to select d max and St satisfying 


T ^max ^max- 


Li((M + D m x )6t H--—) H——b L 2 ( Xv max St + 


Equivalently, by virtue of (7.6) and (7.7), it is required that 
dL 


2 pSt ■ 


-‘'max 

2 


2 st 


+ 2 X^5t + 


2 

d \ r 


<(1 — A )' 


< P~ 


(1 - X)St - ASt 2 
35t +1 

By evaluating the derivative of d max (-) with respect to St we obtain that 

d max = 0 ((1 - A) - 8St){3St + 1) - 3((1 - X)St - 45t 2 ) = 0 

- 12<5f 2 - 85t + (1 - A) = 0 


Hence, we obtain the time 

- = 8 - V4 + 3(1- A) = —2 + a/4 + 3(1 — A) 

—24 ~ 6 

corresponding to the maximum possible diameter 

- (1 - A )St - 4 St 2 

max “ P 3ft+l 

For the simulation results, we select the distance p = 10 and the radius of the circular domain 
R = 10. We also assume that the agents 1, 2, 3 and 4, are initially located at x\q = (5, —3), 
£20 = (5,3), a ; 30 = (0,6) and X 40 = (—4,6), respectively. Thus, it follows that agents 1, 2, 
and 2, 3, satisfy the requirement on their initial relative distance. In the sequel we will focus 
on the behaviour of the system for times t £ [0,2]. Given this time interval and a selection of 
the planning parameter A £ (0,1), we choose the time step St as the largest possible time step 
not exceeding St above, in such a way that the number of time steps NT := ^ is a positive 
integer. We also choose the largest possible cell diameter c? max corresponding to St and consider 
a square grid in R 2 . Each square has side length d , where d is the largest number not exceeding 
-^d max , such that the quotient i s an integer. Thus, we can form a cell decomposition of the 
circular domain D by defining as a cell each square in the grid which has nonempty intersection 
with D. In Figs. 3, 4 we have plotted (half of) the grid lines, in order to illustrate how the grid 
is affected by the choice of A. We next consider two cases for the motion of agent 2, which is 
unaffected by the coupled constraints. 

Case I: It holds v 2 {t) = v 2c , Vf £ [0, 2], with v 2c = (—3, —3). 
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Case II: It holds v 2 {t) = v 2c + v 2 d(t), Vt € [0,2], with v 2c as above and v 2 d G lid, where Ud 
is the set of all piecewise continuous functions v : [0,2] -A- R 2 that satisfy v(t) = 7(i)( 2 > — 2 )> 
with —1 < 7 (t) < 1 for all t G [0, 0.9] and 7 (t) = 0, for all t G (0.9, 2]. 

Notice that in Case I we consider a pre-specified path for agent 2, by selecting a constant 
control, whereas in Case II we allow for the possibility to modify this path and superpose a 
motion perpendicular to it (up to certain bound) over the time interval [0,0.9]. Furthermore, 
in both cases the magnitude of v 2 (-) is bounded by u max (= 5). 

For Case I, we assign reachability goals to agents 1, 3 and 4 which should be fulfilled at 
the end of the time interval [0,2], given the selected path for agent 2. Specifically, we want 
agents 1, 3 and 4 to reach the corresponding boxes in the workspace that are depicted in 
Fig. 3. First, we sample the trajectories of 2 at the time instants kS t, k = 0, 1..., NT and 
specify the sequence • • • l 2 T corresponding to the cells Si K with x 2 (n6t, x 20 ) G Si K . Then, 
we exploit the individual controlled transition systems of agents 1 and 3, in order to determine 
(an underapproximation of) their reachable cells for the given sampled trajectory of agent 
2. In particular, by denoting as l® the index of the cell where the initial state xio of agent 1 
belongs, we can evaluate the indices of its reachable cells at time uSt as Ql = Post^Q” , l 2 _1 ), 
k = 0,1,..., NT, where Q ° := {1°} and we have used the notational convention Postal; l 2 ) := 
Postal; ( l\,l 2 )) (recall that {h,l 2 ) stands for a cell configuration of agent 1 ) and the definition 
PostJ(Qi; Z 2 ) : = U 2 ieQiP 0 S ti(^i; ^ 2 )- The approach followed in this case is possible because agent 
2 is decoupled from the other agents and the individual transition system of agent 1 depends 
only on the cell indices of agent 2. Similarly, we can evaluate the reachable cells of agent 1 and 
check whether it fulfils its reachability task. Next, by computing the reachable cells of agent 3 
which lie in its target box at the final time step NT, we calculate the backward reachable cells 
of the agent in order to encode the discrete trajectories which fulfil its reachability goal. Then, 
we exploit the individual transition system of agent 4 in order to determine its reachable cells 
for all the possible trajectories of agent 3 that satisfy its reachability task. The corresponding 
simulation results are depicted in Fig. 3 for A = 0.2 (left) and A = 0.3 (right). The figure also 
illustrates the effect of the parameter A in the accomplishment of the reachability goals, since 
for A = 0.2 only agent 3 reaches its target box, whereas for A = 0.3 all agents achieve their 
corresponding task. 
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Figure 3. Reachable cells of the agents for (i) A = 0.2 and (ii) A = 0.3. 
Agents 1, 3 and 4 are initially located at the bottom right, top center and top 
left of the illustrated workspace, respectively. The circles denote the sampled 
trajectory of agent 2 as determined by Case I and the boxes the corresponding 
target sets of agents 1, 3 and 4. The union of all discrete paths of agent 3 
which end in its target box are highlighted within the union of its reachable 
cells. 


For Case II, we exploit the individual controlled transition system of agents 1, 3 and 4 in 
order to obtain (an underapproximation of) the cells these agents can reach, irrespectively of 
the choice of v 2 d for the free input of agent 2. In particular, we define the finite cell sequence 
{Q 2 }«;e{o,i,...,AiT} as Q% = {l £ 2 : 3 v 2 ,d £ U d with x 2 (nSt, x 20 -, v 2c + v 2d {-)) £ Si}. Also, we 
inductively define for k = 0,1,..., NT the sets Q} = U iig g^i n Z2g g«-i Postal; l 2 ), Q% = 

^h&Qr 1 nj 2 eQJ -1 p °stl(Z 3 ; h) and Q% = U^qj-i n, 86 Q j-i Postwith Q° = {/?}, 
Q 3 = {Z 3 } and Q% = {^ 4 } (we use the same notational convention as above for the operators 
Post)’('), and the notation l ° for the initial cells of the agents i = 1,3,4). Next, consider any 
selection of sequences l®l\ ■ ■ ■ ^( VT , l®l\ ■ ■ ■ l^ T and l%l\ ■ ■ ■ l ^ T , of agents 1, 3 and 4, that satisfy 
li £ Qi, l}} £ Q 3 and l 1 } £ Q f, respectively. Then, by taking into account the definition 
of the sets Q%, i = 1,3,4, the definition of the individual transition systems of agents 1, 3, 
4, and the particular coupling between the agents in this example, we arrive at the following 
conclusion. For each agent 1, 3 and 4, it is possible to assign a sequence of control laws, 
such that each corresponding agent will reach the cells with indices 1 4 , 1 3 and 11 at time uSt, 
respectively, for any selection of the input v 2 d of agent 2. In Fig. 4 we illustrate the union 
of the reachable cells of agents 1, 3 and 4 for A = 0.3 and A = 0.4, respectively. Notice that 
the underapproximation of agents’ 1 and 3 reachable cells increases with the selection of the 
larger parameter A, namely, with the exploitation of a larger part of the free input for planning. 
However, the same observation does not necessarily hold for the reachable cells of agent 4. The 
reason why the area covered by the reachable cells of agent 4 remains approximately the same, 
is that the corresponding area increases for agent 3 for larger values of A. Thus, although 
the reachability properties of agent 4 are improved, this is compensated by the fact that each 
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illustrated transition of agent 4 to a certain cell needs to be possible for an increasing number 
of different positions of agent 3. 



Figure 4. Reachable cells of the agents for (i) A = 0.3 and (ii) A = 0.4. 
Agents 1, 3 and 4 are initially located at the bottom right, top center and top 
left of the illustrated workspace, respectively. The circles denote the nominal 
sampled trajectory of agent 2 and their nearby cells represent the cells where 
agent 2 can lie at the sampling times, for all possible inputs of Case II. 


The code for the simulation results has been implemented in MATLAB and the worst case 
running time for the illustrated results is of the order of 45 minutes, on a PC with an Intel(R) 
Core(TM) i7-4600U CPU @ 2.10GHz processor. 


8 . Conclusions 

We have provided a decentralized abstraction methodology for multi-agent systems and quan¬ 
tified bounds on the space and time discretization in order to obtain meaningful transition 
systems with multiple transition possibilities. The abstraction framework is based on the de¬ 
sign of hybrid feedback control laws that take into account the agents’ coupled constraints and 
guarantee the implementation of the discrete transitions by the continuous time controllers. 

Ongoing work includes the improvement of the acceptable choices of d max and St in order 
to obtain coarser abstractions and reduce the size of each agent’s transition system. Another 
possible direction for complexity reduction is the modification (localization) of the current 
framework through an event based online abstraction with updated choices of d max and St. 
Finally, it should be noted that while this report provides informal indicators of how the results 
can be used for planning, we are currently formalizing a distributed planning methodology from 
high level specifications that builds on the derived abstractions. 
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9 . Appendix 


In the Appendix we provide the proofs of Propositions 4.8 and 4.9 and certain technical parts 
from the proof of Proposition |6.2| 


Proof of Proposition Indeed, consider a final cell configuration Y = (l\ ) • ■ • i I'n) as in the 
statement of the proposition. By the definitions of the operators Postj(-), i £ Af and the 
transition relation of each corresponding agent’s individual tra nsiti on system, it follows that 

each transition l, p -±Yf j s we n posed in the sense of Definition 


P r i(f) 


4.3 


Thus, given the initial 


cell configuration pr, ; (1) and the well posed transition li — A l[ for each agent, it follows from 
Definition 4.3 i) that we can pick for each i £ Af a control law /Ci >pr (I) (•) that satisfies Property 
(P) and a vector Wi £ W, such that Condition (C) holds. Next, notice that for each agent i the 
projection of the initial cell configuration 1 = (Zi,..., Zjv) is pr, ; (l), namely, the cell configuration 
according to which the feedback law fcj lPr .(i)(-) was selected, and that the feedback law satisfies 
Property (P). Thus, it follows from Condition (C) that the solution of the closed loop system 
is well defined on [0, <5t], and that for each i £ Af, the i-th component of the solution satisfies 

dO. □ 


Proof of Proposition \4-9\ Proof of (i). Let 'w t £ W , i £ Af and x(0) £ D N with ajj(0) £ Si it 
i £ Af be the initial condition of the closed-loop system. Then, it follows from the local 
Lipschitz property on the functions /)(■) and the corresponding property on the mappings 
^s.prji) (•; x iOi w i) provided by (PI), that the dynamics of the closed loop system are given by 
a locally Lipschitz function on D N . Hence, there exists a unique solution x(-) = x(-,a;(0)) to 
the initial value problem, which is defined and remains in D N for all times in its right maximal 
interval of existence [0,T max ). We proceed by proving that each component a;*(•), i £ Af of the 
solution satisfies 

Xi(t) £ ( S U +B(R max)) ^ D, Vi £ [0, min{T max , T}). (9.1) 

Indeed, suppose on the contrary that ( |9.1| ) is violated, and hence, by taking into account that 
Xi(t) £ D for all t £ [0, T max ), that there exists i £ Af and a time T with 

T £ (0, min{T max , T}) and x L (T) £ S L + B(R max ). (9.2) 

By recalling that Xi(0) £ Si t , i £ Af, we may define 

r := max{t £ [0,T] : a;*(*) £ c l(S h + B(R max )), 

Vt£[0,t\,i£Af}. (9.3) 


Then, it follows from (4.9), (9.2) and ( |9.3| ) that there exists £ £ Af such that 

xi[r) £ d(Si e + B(R max )) 


(9.4) 


and that 

t <T <T < St. 


It also follows from (9.3), (|4.9[), (|9.5[) and Property (P2) that for all t £ [0,r] it holds 


(9.5) 


l^,pr^(l) (i; Xg(t'), (t) , XgQ, W(f) \ — ^max 


( 9 . 6 ) 
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Hence, we get from (|3.l|), (3.2), (4.1), (|4.8|), (9.6) and (9.5), which implies that r < St, that 


\xe(T)-xeo |< / \fe(x e (s),x jW (s))\ 

Jo 

+ | h, pr^(l)(s, xe(s),-Xj(t)(s)\ X£o, wg)\ds 
< [ (M + v max )ds < St(M + u max ) = Rn 


(9.7) 


In order to finish the proof of (9.11 we exploit the following elementary fact. 

Fact I. Consider a nonempty set S C and a constant R > 0. Then for every x £ d(S+B(R)) 
it holds \x — y\> R,\/y £ S. 

Proof of Fact I. Indeed, suppose on the contrary that there exists y £ S with \x — y\ < R — e 
for certain s > 0. Then for all x £ int (B(x;e)) we have 

\x — y\ < \x — a:| + \x — y\ < e + R — e = R, 

and hence, x £ S + B(R) for all x £ int(H(a;;e)), which implies that x £ d(S + B(R)) and 
contradicts our statement. < 


By exploiting Fact I with S = Si t , R = R m ax, y = Xfo and x = xi{t) we deduce from (9.7) 
that xi{t) £ d(Si t + B(R m ax )) which contradicts (9.4), and provides validity of (9.1). 

We now prove the following claim: 

Claim I. It holds T max > T. 

>Proof of Claim I. Indeed, suppose on the contrary that 

Pmax < T. (9.8) 

For each i £ J\f let Vi : [0, oo) -> R n be a piecewise continuous function satisfying 

Vi(t) = k i}PTi (i)(t,Xi(t),-x.j(t);Xio,Wi),yt £ [0,T max ). (9.9) 


Notice that due to H.9P and (9.8) we have that T, 


< minjdt, min{T(a; i o, w?:) : i £ A/”}}, and 
: , Vt £ [0, T max ). Hence, we may select 

). Thus, if 

we denote by £(•) the solution of (3.1 )-([3(2| with free inputs Vi(-), i £ N and the same initial 


thus, we get from ( |9.1[ ) and (P2) that |ui(t)| < v n 
Vi(-) to satisfy \i’i(t)\ < u max , Vf > 0 (select for instance Vi(t) = 0 for t > T, 


condition with x(-), it follows from the Invariance Assumption (IA) that £(f) is defined and 
remains in D N for all t > 0. Furthermore, it follows from standard arguments from the theory 
of ODEs that £(f) = x(t),Vf £ [0,T max ). Hence, since ^(f) belongs to a compact subset of D N 
for all t £ [0,T max ], the same holds for x(t) on [0,T max ). The latter contradicts maximality 
of [0,T max ) since by (9.8) and (4.9) it holds T max < T < min{T(a,’io, Wi) '■ i € A/”} and the 


mappings fcj iPr .(i)(-) are defined for t £ [0, min{T(a:io, Wi) '■ i £ A/”}). Hence, we have shown 
Claim I. < 

From Claim I, it follows that x(t) is defined and remains in D N for all t £ [0, T) and that 


(9.1) holds for all t £ [0, T). Thus, by applying the same arguments with those in the proof 
of Claim I, we can determine a continuous function £(•) with £(t) = x(t) for all t £ [0, T) and 


£(T) £ D n , which establishes (4.10). 


Proof of (iia). In the case where (P3) also holds, and hence by ( |4.9| ) we have that T = St, it 
follows from part (i) of the proposition and standard arguments, that the solution x{-) is defined 
on [0, T max ), with T max > St. From the latter, we conclude that x{t) £ D N for all t £ [0, St]. 
Moreover, since T max > St = T, it follows that (9.1) is satisfied for all t £ [0, St). The latter, by 
virtue of (P2), (P3) and continuity of x(-) implies (4.11). 
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Proof of (iib). By exploiting the result of part (iia) of the proposition and defining Vi(t) = 
fcj pr (i) ( t , Xm, Wi), Vf € [0, St) we can extend Vi(-) to a piecewise continuous function 

on [0, oo) which satisfies (3.4). Hence, by applying the same arguments with those in the proof 
of Claim I, we conclude that the solutions x(-) of (3.1 )-<3.2>, (4.8) and £(•) of (3.1)-(3.2) (with 
input v(-)) coincide on [0, <5t]. □ 


Proof of the fact that the intervals provided by Cases I, II, III in Proposition |6.2| 
are well defined. 


The fact that (6.6) and (6.7) are well defined is straightforward. We proceed by defining 


S'(m) : = 


2(A(1 - A)/x - 2A 2 )< ax 
g 2 ML 


(9.10) 


Then it follows that 


g(/x) > 0 


A(1 - X)n - 2A 2 > 0 


A* > 


2A 

1 - A' 


(9.11) 


Furthermore, 


^ , d f { 1 - X)n~2X 

sgn(sw) = sgn ( — 


= sgn 


p 

(1 — A)/i 2 — ((1 — X)fi — 2X)2fi 
M 3 


= sgn((l - X)n - 2((1 - X)n - 2A)) = sgn(-(l - A )/i + 4A), 


which implies that 


4A 


sgn (g{/j))[> or =]0 n[< or 


(9.12) 


From (9.11) and (9.12) we get that 


0< JW < s(j i ^)= I 4; 


4AV 2 


A 2 

(1-A) 


-ML 


(1-A) 2 <, 

4 ML 




2A 4A 


1 — A’ 1 — A / V 1 — A’ 


U 


4A 


(9.13) 


which implies that (6.8) and (6.10) are well defined. Also, by considering the function 


h (d max ) • — 


(1 - A)U max + \/(l ~ A) 2 ug, ax - 4 MLd n 
2 ML 


(9.14) 


we deduce that 


ix < 


4A ii 

1 — A 2Au m 


< 


(1 A)u max „ 


2 ML 


< h (d max ),Vd max e 0, 


(! A ) 2 <ax ^ (9.15) 


4ML 
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and thus, that (6.91 in Case II is well defined. Finally, in order to also show the latter for the 
interval (6.9) in Case III, we note that due to (9.10) and (9.14) it holds that 


h"(gM) = 


(1 - A)W + y/(l- A) 2 < ax - 4 ML 2{m1 


2 ML 


(1 - A)u max + yUax ((1 - A) 2 - 8A(1 - A) + 16 (A)') 


2 ML 


(1 A)u max 4~' 




(1 — A) + (l-A)-f 


2 ML 

Hence, by taking into account that 


2 ML 


(1 — A) — — > 0 
I 4 


M > 


4A 

1 — A 


and that in Case III it holds fi > 44 we get that 

2 ((! - A) + Wmax _ ((i _ A)/x - 2A)r 


h"(3(n)) = 

Furthermore, we have that 


2 ML 




-.9 (A 4 ) = 


M 2(A(1 - A )/j - 2A 2 X. 


UML 

((1 A)/x 2A)u max 


2Ai.’ ma x 2Aw max fi 2 ML jiM L 

Thus, from (9.13), (9.17), ( |9. 18 ) and the fact that h"{■) is decreasing, it follows that 


h"(dmax) > 


2 \v n 


^maxj \/d 

max e (0,g(/i)], 


(9.16) 


(9.17) 


(9.18) 


(9.19) 


which in conjunction with (9.10), implies that (6.9) is well defined. 


Proof of the fact that for all Cases I, II and III (6.18)-(6.21) are satisfied in the proof 
of Proposition |6.2| . 

Case I: 0 < /i < 44 
By defining 


h (drnax) • 


we obtain that 


(1 - A)U max - \/(l - A) 2 Um ax - 4 MLd n 
2 ML 

1 


(9.20) 


^ (^max) 


Hence, 


x /(l-A) 2 < ax -4MLd max ' 
h!{•) is positive and strictly increasing for 0 < d max < 


(l~A) 2 r4 

AML 


(9.21) 


and furthermore 


h'( 0 ) - 


1 


(1 - A)v n 


h'(0) = 0 . 
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The latter in conjunction with (9.21) implies that 


h (t^max) ^ 


1 


> 


1 


(1 - X)v max “““ M + v n 
Also, we have that for 0 < fj, < ^4- it holds 


dmax 7 ^ dmax F: I 0 5 


(1 — X) 2 V . 


2,,2 

max 


4 ML 


2Xv n 


< 


2A 
1 —A 


2Ai; n 


" dniax — 


1 


(1 - X)v D 


(9.22) 


(9.23) 


Thus, it follows from (6.6), (6.7), (9.20), (9.23) and ( 9.22| ) that (6.18)-(6.21) are fulfilled (see 
also Fig. 5). 


St 

4 



_2_ 

(I-A)jj 


(1 A)u ma , 

_ H _ ,] 

2A» max 

_ 1 _ c 

AT-t-'U m ax 

dma.x 


FIGURE 5. Case I. Feasible d max — St region for 0 < fi < 


Case 

From 


II: 


— < u < — 

1 —A ^ ^ 1 - 


1 —A ‘ 


(9.10) and ( 9.20| ) we obtain that 




(1 - A)t! max - yj{ 1 - A) 2 < ax - 4Af£— (1 X ft M 2 L X ~ )v; ^ 

2 ML 

(1 - A)w - W< ax ((1 - A)^ - 8^(1 - A) + 16 (£) 2 ) 


2 ML 


(1 A)'Umax %ax^ 

31 =■ 

i—1 

2 

’(I - A) - 

(1-A) 

2 ML 

2 ML 


4A 


^max 


Hence, by taking into account (9.16) and that y^ < F < y^, we get that 

h ,,- ( u _ 2 ( (1 ~ A) ~^)^ max _ ((1 - X)n - 2Xn)v max 
K9W) 2 ML /iML 


(9.24) 


2 ML 
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Thus, by the fact that < fi < it follows from ( 9.13 ), ( 9.18 ), ( 9.21 ) and ( 9.24 ) that 


h (d max ) 5: ^ a ^max^Vdmax G (0,^(/i)] 


(9.25) 


and 


d (f^max) ^ o \ ^ max 5 Vd 

max e (s(d) 

2A'U max 


(l-A) 2 ^i 


2„,2 

max 


'1M L 


(9.26) 


It also holds that 


^ 1 A = ' > dmax ^ /. » \ d max > dmaxj Vd max A 0. (9.27) 


1 — A 2Az\, 


(1 - A)u n 


M + v n 


Hence, it follows from (6.8[)-( |6.7[ ), ( |9.10[ ), ( |9.25[ ), ( |9.26| ) and ( |9.27[ ) that ( |6.18[ )-( |6.21[ ) are fulfilled 
as well (see also Fig. 6). 


St 

A 



FIGURE 6. Case II. Feasible d max — St region for < fr < yzj; 


Case III /r > 


4 A ^ — If— A - ) " m " we have 


In this case, notice that for /z = and d max = 4ML 
and thus, we deduce from (9.21) that 


2Ad„ 


dmax — d (dmax)j 


d > 


4A 

1 - A 


2Az; n 


> h'(d, 


Vd n 


max ) •> v u max 


e o 


(l-A)X, 

4 ML 


(9.28) 


Hence, it follows from (|6.8|), (|6.9|), (|9.10|), (|9.13|) and (|9.28|) that (|6.18|)-(|6.21|) hold (see also Fig. 

7). 
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St 



FIGURE 7. Case III. Feasible d max — St region for p > 


10. Acknowledgements 

This work was supported by the H2020 ERC Starting Grant BUCOPHSYS and the Swedish 
Research Council (VR). 


References 

[1] A. Abate, A. Tiwari, ans S. Sastry, Box invariance in biologically-inspired dynamical systems, Automatica 
45(7):1601-1610, 2009. 

[2] D. Angeli and E. D. Sontag, “Forward completeness, unbounded observability and their Lyapunov charac¬ 
terizations,” Syst. and Control Lett., vol. 38, pp. 209-217, 1999. 

[3] E. Asarin, T. Dang, and A. Girard, Hybridization methods for the analysis of nonlinear systems, Acta 
Informatica, 43(7):451-476, 2007. 

[4] C. Baier and J. P. Katoen, Principles of Model Checking, Cmbridge, MA: The MIT Press, 2008. 

[5] D. Boskos and D. V. Dimarogonas, Decentralized Abstractions for Feedback Interconnected Multi-Agent 
Systems, in Proceedings of the 54th IEEE Conference on Decision and Control, 282-287, 2015. 

[6] D. Boskos and D. V. Dimarogonas, Robust Connectivity Analysis for Multi-Agent Systems, in Proceedings 
of the 54th IEEE Conference on Decision and Control, 6767-6772, 2015. 

[7] M. E. Brouke and M. Gannes, Reach Control on Simplices by Piecewise Affine Feedback, SIAM Journal on 
Control and Optimization, 52(5):3261-3286, 2014. 

[8] P.E. Caines and Y.J. Wei, The hierarchical lattices of a finite machine, Systems and Control Letters, 
25:257-263, 1995. 

[9] Y. Chen, X. C. Ding, A. Stefanescu, and C. Belta. Formal approach to the deployment of distributed robotic 
teams, IEEE Transactions on Robotics, 28(1):158—171, 2012. 

[10] G. E. Fainekos, A. Girard, H. Kress-Gazit, and G. Pappas, Temporal Logic Motion Planning for Dynamic 
Robots, Automatica 45(2):343-352, 2009. 

[11] S. G. Loizou and K. J. Kyriakopoulos, Automatic Synthesis of Multi-Agent Motion Tasks Based on LTL 
Specifications, in Proceedings of the 43rd IEEE Conference on Decision and Control, 153-158, 2004. 

[12] A. Girard and S. Martin, Motion Planning for Nonlinear Systems using Hybridizations and Robust Con¬ 
trollers on Simplices, in Proceedings of the 47th IEEE Conference on Decision and Control, 239-244, 2008. 

[13] A. Girard and S. Martin, Synthesis for constrained nonlinear systems using hybridization and robust con¬ 
trollers on simplice, IEEE Transactions on Automatic Control, 57(4):1046-1051, 2012. 

[14] A. Girard, G. Pola, and P. Tabuada, Approximately bisimilar symbolic models for incrementally stable 
switched systems, IEEE Transactions on Automatic Control, 55(1):116-126, 2010. 















40 


D. BOSKOS AND D. V. DIMAROGONAS 


[15] E. A. Gol, X. Ding, M. Lazar, and C. Belta, Finite Bisimulations for Switched Linear Systems , IEEE 
Transactions on Automatic Control, 59(12) :3122-3134, 2014. 

[16] L. Grime, Asymptotic Behavior of Dynamical and Control Systems under Perturbation and Discretization, 
Berlin, Germany: Springer-Verlag, 2002, vol. 1783. 

[17] L. Habets, P. Collins, and J. van Schuppen, Reachability and control synthesis for piecewise-affine hybrid 
systems on simplices, IEEE Transactions on Automatic Control 51(6):938-948, 2006. 

[18] L. Habets and J. van Schuppen, Control of piecewise-linear hybrid systems on simplices and rectangles, 
Lecture Notes in Computer Science: vol. 2034 Hybrid Systems Computation and Control, 261-274, Springer, 
2001 . 

[19] M.K. Helwa and P.E. Caines, In-Block Controllability of Affine Systems on Polytopes , in Proceedings of 
the 53rd IEEE Conference on Decision and Control, 3936-3942, 2014. 

[20] I. Karafyllis and J.-P. Jiang, Stability and Stabilization of Nonlinear Systems , London, U.K.: Springer- 
Verlag, 2011. 

[21] M. Kloetzer and C. Belta, A Fully Automated Framework for Control of Linear Systems from Temporal 
Logic Specifications , IEEE Transactions on Automatic Control 53(l):287-297, 2008. 

[22] G. Pappas, Bisimilar linear systems , Automatica, 39(12):2035-2047, 2003. 

[23] G. Pola, A. Girard, and P. Tabuada, Approximately bisimilar symbolic models for nonlinear control systems, 
Automatica, 44(10):2508-2516, 2008. 

[24] G. Pola, P. Pepe, and M. D. di Benedetto, Symbolic Models for Networks of Discrete-Time Nonlinear 
Control Systems, in Proceedings of the American Control Conference, 1787-1792, 2014. 

[25] G. Reissig, Computing Abstractions of Nonlinear Systems, IEEE Transactions on Automatic Control, 
56(ll):2583-2598, 2011. 

[26] Sontag E, Mathematical Control Theory, 2nd ed. New York: Springer-Verlag, 1998. 

[27] P. Tabuada, Verification and Control of Hybrid Systems, A Symbolic Approach, New York: Springer, 2009. 

[28] Y. Tazaki and J. ichi Imura, Bisimilar finite abstractions of interconnected systems, in Hybrid Sys¬ 
tems: Computation and Control, 4981:514-527, 2008. 

[29] A. Tiwari, Abstractions for hybrid systems, Formal Methods in Systems Design, 32(l):57-83, 2008. 

[30] A. Van der Schaft, Equivalence of Dynamical Systems by Bisimulation, IEEE Transactions on Automatic 
Control, 49(12):2160-2172, 2004. 

[31] M. Zamani, M. Mazo, and A. Abate, Symbolic Models for Networked Control Systems, arXiv: 1401.6396, 
2014. 

[32] M. Zamani, G. Pola, M. Mazo, and P. Tabuada, Symbolic Models for Nonlinear Control Systems Without 
Stability Assumptions , IEEE Transactions on Automatic Control, 57(7): 1804-1809, 2012. 

Department of Automatic Control, School of Electrical Engineering, KTH Royal Institute of 

Technology, Osquldas vag 10, 10044, Stockholm, Sweden 

E-mail address: boskos@kth.se 

Department of Automatic Control, School of Electrical Engineering, KTH Royal Institute of 

Technology, Osquldas vag 10, 10044, Stockholm, Sweden 

E-mail address: dimos@kth.se 



